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Third Paper on Tensor Analysis.” 
By G. Y. 


This paper deals with some questions of the integral calculus of tensor 
fields using the method and notions introduced in the foregoing papers. 

Before we state, by way of introduction, what questions are treated here, 
it will be well to point out a question which will not be discussed in the 
body of the paper. The differentiation of a tensor field of rank & leads to a 
tensor field of rank & +1 (First paper §§4 and 9; second paper § 18). 
It.is natural, given a tensor field of rank & +1, to ask whether it may be 
considered as a differential of a field of rank & and, if so, to try to find this 
tensor field. In the case k —0 the given field would be a field of rank one, 
i, e. a vector field, and the questions whether this field is a gradient and, 
if so, how to find the scalar field of which it is the gradient, are. answered 
in a way entirely analogous to that followed in a Euclidean space. If; how- 
ever, k > 0, the situation becomes radically different: let us consider. the 
case k = 1, for instance, using the customary coordinate notations; if ¢; is a 
tensor field of rank one, if ¢;,; is its first and ¢i,j, its second differential, 
it is well known that 


i, ik — = Ri. jx: hp, 


where £‘;. jx, is the Riemann tensor of the space. If therefore a tensor field 
Wi; of rank two is given we see that, if it is the differential of a tensor field 
¢; it must. satisfy the relation 


Wijk — = hp; 


which in a general space determines the primitive field completely. We have 
thus a true “ integration by differentiation ” and the primitive field is deter- 
mined uniquely, there are no constants of integration. We shall not discuss 
in this paper these interesting questions, which, it seems, are also import- 
tant from the physical point of view; we shall not deal here with the research 
of primitive functions, not with questions which correspond to the study of the 
indefinite integral but with the process of integration, with the generalization 
of the definite integral. : 

We shall follow the system adopted in the first two papers, beginning 


* The two first. papers appeared in this Journal, Vol. XLVI, No. 2, April,. 1924, 
pp. 71-94 and Vol. XLVII,. No. 1, January, 1925, pp. 11-24. 
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with the study of ‘the simpler case of ‘a surface -in the ordinary space; the 
indicator, as it ought to be, shows the way; it shows that integration leads 
away from the surface into the space; we then take up the general case of q 
tensor field on the surface, show what integration becomes in this case and 
deduce the integrability conditions or the conditions for the independence of 
the integral of the path; besides a differential condition which is the exact 
counterpart of the usual condition we have here additional algebraic condi. 
tions. Before entering upon the discussion of the generalization of these 
things for higher spaces we show how the Riemann tensor can be introduced 
in a new way which in many respects resembles the introduction of the 
indicator in the’ case of a surface. The idea of tensor area seems to prove 
useful in this connection as well as elsewhere. After we introduce surface 
integration and volume integration for higher spaces and discuss the integra- 
bility conditions following the way adopted in the case of the surface we treat 
the Riemann tensor and the contracted Riemann tensor as examples. In the 
last section we indicate a method of obtaining integral expressions without 
leaving the curved space.* 


§ 21. Line Integrals on a Surface. 


We remember that the (negative) indicator was obtained in our first 
paper (7.5) as 


—h, 


b—a 
lim with 
where b and a were unit normal vectors at the points B and A respectively. 
We may consider the totality of unit normal vectors as a vector field (which, 
however, does not belong to the surface). If we introduce the notation 
b=f(B) we have 


so that the negative indicator appears as the differential of this vector field 
(compare 4.1). We may also say that the indicator gives the infinitesimal 
change of the unit normal vector which corresponds to the infinitesimal dis- 
placement h on the surface. If we write dA for h, we have 


da = —s(dA). 


* The contents of this paper was briefly outlined in papers read before the Ameri- 
can Mathematical Society on December 28, 1924 and February 27, 1925. 
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3 the It is natural to consider the finite change of the unit normal vector g—p if 
leads as the integral of —s(dA), 
of a @ | 
a This suggests going outside of the surface if we want to consider 
ndi. integrals. | | 
al Following this suggestion we shall find an operation of integration of il 
which this transition from s(dA) to qg—p is a special case. One way to il 
do this is to consider the integral as the limit of a sum. It is, however, not HI 
b necessary to introduce this process explicitly, as it is possible to reduce the i 
face case which we have to consider to things already known. The hint which we i 
have, viz., that we have to take into account the outside space, will help us ‘ 
sd to make this reduction. 
the The situation is this. If we have a vector field, i. e. a vector v4 at each 
point A, we form the integral Hl 
out i i 
fva.dA 
where v4.d@A is the scalar product of the vectors v4 and dA; this can be con- | | 
sidered as the limit of the sum 
rst v4.(B— A) + ve.(C—B) + 10.(D—C) +: 
where A, B, C, D,- ~~ are points along a curve, and in coordinates has the 
familiar form 
+ vydy + vzdz. 
y: Instead of considering a vector v we may consider (property 1, § 3, First | 
h, Paper) the scalar linear function—tensor of rank one—¢(x) =v.2; our 
” integral takes then the form 
f¢(dA). 
If we have, say, a tensor of rank two, ¢(z,y) and our surface happens to be 
1 a plane, we may give to the second argument y a fixed value and form, without 
] any difficulty, the integral 
So(daA,y), 
(for every given path of integration) a function of y, to be precise a tensor 
of rank one whose argument is y, instead of a number which we obtained 
before—and the same can be applied to tensors of higher ranks. But, if we 
have a curved surface, we encounter the same difficulty as in the case of 
differentiation (§ 9): we cannot give to y a fixed value. We therefore resort 
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to the following procedure which is of fundamental importance for all thy 
follows: 

We take a fixed, but arbitrary vector c of the space which contains oy 
surface; we consider in every point A of the surface a vector ca which is th 
projection of c on that bundle (on the tangent plane at A) and in formiy 
the integral, we give to y at each point A the value Ca, viz., 


f $(dA, Ca). 


We obtain thus (for every path of integration) a number, say y. No 
we make c variable (this corresponds to the restoration of variability toy 
in the case of the plane). The number y becomes a function of the vector, 
and, as we easily see a linear function (using the fact that the projection 
the sum of two vectors is the sum of their projections and the definition of 
a linear function in § 3, First paper). We thus obtain as a result of integn. 
tion of a tensor field ¢(z, y) of rank two (for a given path) a tensor of rank® 
one but this tensor is a tensor of the three dimensional containing space 
We may, again, consider instead of the tensor y(c) the vector w defined by the 
relation y(c) —w.c (§3) and instead of the tensor of the second rank 
(x,y) we may consider, at each point of the surface, the equivalent linea 
vector function g(x) which is given by g(x).y—¢(2,y) and say that the 
result of integration of a linear vector function, which belongs to a surface, 
along a path on this surface is a vector of the containing space. 

The indicator being a linear vector function, and the difference between 
the unit normal vectors at the ends of a path being a vector of the containing 
space, we see easily that we have actually found the generalization which we sé 
out to find. 

The further generalization for fields of higher ranks is obvious. In order 
not to complicate the writing we shall only mention the case of rank three; 
we take in this case two arbitrary vectors of the containing space ¢ and 4, 
consider their projections c4 and qa, form as before the tensor field of rank 
one $(dA, ca, ga), integrate this along a line T and obtain a tensor of the 
containing space y (c, q) which has to be considered as the integral: 


If the given tensor field is not symmetric we may also form two other 
integrals 


f ¢(c,dA,q) and ff 
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§ 22. Integrability conditions. 


ry : In the case of the indicator we know, from the geometrical meaning of 
form the integral f s(dA), which gives the change of the unit normal vector along | 


the curve I, that if two curves [, and I, have the same end points Wi 


s(dA) = f s(dA), | 


Ti Ts 


: provided the curves may be continuously transformed one into the other with- a 
ty out leaving the surface, because each of the two integrals represents the dif- i 
he ference between the normal unit vectors at the ends of the path. This sug- i 
of i} 


gests considering the general question of when, given an integral of a tensor 
field of rank two, it depends only on the end points of the path of integration 
and not on the path itself or when the integral taken over a closed reducible 


contour vanishes. In order that this be so f ¢(da,c4) must be independent i 
of the path for every value of c.’ If we put, for a given c, $(2, ca) = (2) 


we know that the condition of independence is 


22.1 w(x, h) =’ (h, | 


where w’(z,h) is the differential of the tensor field w(x). (It would lead 
us too far if we desired to give here the vector deduction of this result 
which is well known in coordinate form). We have, therefore, to find the 
differential of the tenSor field $(x4,c4) (2a), where ca is the projection 
of a constant vector of the containing space. We know (First paper, § 9) 
that the differential of w(x) or, to be more explicit, w(A,2z4) is given by 


o(B, 4B) w(A, La) 
r 


with as — ha 


w’(A, ha) == lim 


where 24.z is the projection of z4 on the bundle B. We therefore have 


22.2 w’(A LA ha) = lim $(B, Cr) — ¢(A, LA; Ca) 
» VA; 


This is not the differential of the tensor field of the second rank 
$(A, 24, ya) for ya ca. To have the differential we must have in the first 
term of the numerator c,.z, the projection of c4 on the bundle B, instead of 
cg, the projection of c on the bundle B. However, this coincides with the 


other 
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differential $’(A, 2a, ca), if ¢ coincides with ca, i. e. if the vector belongs 
(or is parallel) to the bundle A. In this case 22.1 takes the form 


22.3 CA; ha) (ha, CA; La) 


Since 22.1 must hold for all values of c, 22.3 must be satisfied, it is q 
necessary condition of integrability; in order to find sufficient conditions we 
will have to consider the case when c does not belong to the bundle A. Before 
we do this we notice that, applied to the indicator, equation 22.3 gives us the 
Codazzi equations of which we have thus found a derivation which seems to 
be novel. 

Coming now to the general case we may present c as the sum of a vector 
ya of the bundle A and of a vector which is normal to this bundle and which 
we may write pa, a being the unit normal vector and p a number. We have 
thus 

c=yatpa, Ca=—Ys, Ca= pas; 


22.2 takes the form 

$(B, — (A, ya) + pb(A, ap) 

== 24, Ya; ha) + po(A, La, lim ag/d) 


We have to find now lim ag/d. The vector ag is the projection of the unit 
normal vector at A on the tangent plane at B and is easily shown to be 
a—b(ab) using the reasoning employed in proving formula 9.2. We have 
now 


22.21 w’(A, VA; ha) = lim 


a—b(ab) a—b 


lim = lim lim 


(compare the similar computation in § 12). We have thus 
22.4 w’(A, TA, ha) = (A, LA, YA; ha) po(A, LAs» $(ha) 


Since 22.1 must hold for an arbitrary c we have to consider here y4 and p 
as arbitrary. Making p—0 we obtain once more condition 22.3; making 
ya =0, p=1 we obtain from 22.1 


22.5 $(x, 8(h)) = $(h, s(z)) 


where the letter A is dropped as argument and as index. This together with 
22.3 which we write more conveniently as 
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22.6 ¢’ (2, h) ¢’(h, x), 


gives a set of necessary and sufficient conditions of integrability. 
If we have the integrand in the form of a vector function f(z), the con- 
dition 22.5 becomes 


22.7 f(x) .s(h) =f(h).s(z) 
which may also be written as 
22.71 .h=sf(h).x 


and, in the case when f is symmetric, in the even simpler form 
22.72 sf(x) =fs(x) or sf—fs 


which shows that the function f has to be commutative with the indicator s; 
condition 22.6 takes here the form 


22.8 h) = f’(h, 2). 


In the case of the indicator condition 22.7 is clearly satisfied, both sides be- 
coming identically equal. 

It may not be superfluous to insist on the fact that we have here, besides 
the usual differential condition of integrability, an additional algebraic con- 
dition involving the indicator. If both these conditions are satisfied for a 
tensor field of rank two which belongs to the surface there exists an oo of 
vector fields which are defined at all points of the surface but do not belong 
to the surface, because their vectors instead of being tangent stick out of it. 
All these fields may be obtained from one of them by the addition to all vec- 
tors of an arbitrary constant (space) vector. 

If the integrand is a field of rank three $¢(2,y,z) we have to put 
c=ya + pa, g = 24 + 00; we obtain instead of 22.21 and 22.4 


(A, TA, ha) lim $(B, YA.B; Z4.B) $(A, VA; YA, Za) 


po(A, ap, 24) La, Ya, Op) + pod(A, La, Op, ap) 


= $’(A, ya, Za, ha) + pb(A, 8(ha), 2a) + 06(A, La, ya, 8(ha)) 


+ lim 


the term containing ag twice vanishing in the limit. 
The integrability condition consists of three relations, one differential : 


22.61 y, 2,h) = $’(h, y, 2, x) 
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and two algebraic: 

22.51 s(h), z) = (2), 2) 
and 

22.52 $(x,y,s(h)) = $(h, y, 8(2)) 


§ 23. Line Integrals in Higher Spaces. 


If, instead of a surface of the ordinary space, we consider an n-S. im- 
mersed in an N.-E.S., line integrals of tensor fields can be introduced in exactly 
the same way. In deducing the integrability condition a change has to be 
made, however. The vector c which belongs to the containing space must he 
presented as y4 +a, where a is an arbitrary normal vector; we cannot put 
it in the form pa there being no definite unit normal vector when r >1, 
Formula 22.21 takes the form 


23.1 w’(A, 24, ha) = $’(A, Ya, ha) + (A, Xa, lim ap/d) 


where a is an arbitrary normal vector. We have now to find lim ag/d in 
this case. 

We shall use this occasion and take up again the question of introducing 
the indicator dealt with in §17 of the second paper and then in the added 
remark at the end of it. We consider at each point B of the neighborhood 
of A, B approaching A in such a way that (B—A)/A—ka, a tangential 
vector hg and a normal vector b, the values of these vectors for A being ha 
and a respectively. We have then 


23.2 hg.b =0 ha.a=0 


We suppose that and approach finite limits as A approaches 


zero so that 


This gives 


lim = lim 


Tim 


or in view of 23.2 


23.3 — lim 


Since the left hand side of the equality does not contain hz and the right 
hand side does not contain 6 we conclude that the common value of both sides 
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does not depend on either. This common value is, therefore, only a function 
of the vectors ha, a and, of course, of k4. Let us denote this function which, 
obviously, is a linear function of each of the three vectors, by y(ha, ks,a), so 


that we have 


lim 


23.4 W(ha, ka, a) = — lim 


Using property 3 of §3 (First paper) we may write y also in one of the 
following forms 


23.5 W(ha, ka, @) = ha F(a, ks) =a.8(ha, ka) 


where F(a, ka) is the tangential component of lim ng and s(ha,ka) the 


h 
normal component of aaa oe . The functions § and s have been intro- 


duced in the second paper ( ¥ with the opposite sign). Here F (a, ka) is a 
linear vector function one of whose arguments is a tangential and the other 
a normal vector, while its value is a tangential vector; both arguments of s 
are tangential vectors and its value is normal; y is a scalar trilinear function 
with two tangential and one normal argument. The functions y, F, s are 
only different forms of the same thing. We shall refer to all three of them 
indifferently as to the indicatrix (reserving the name indicator for the function 
$(x) which exists only in the case N—n=1). 

We return now to our question: What is lim ag/A? It is clear, first 
of all, that it is a tangential vector. Denoting by a» the normal component 
of a at the point B we may write 


ap = ad 
and 
a— ap ay. h 


A 


= — lim 


lim 
Making now b = a in the formulas 23.4 and 23.5 we may write 


ha. F(a, ka) = lim — 


and, since lim ag/A is tangential, as remarked, 
23.6 lim ap/A = — F(a, ka). 


There would be some advantage in taking this as the definition of ¥. Rela- 
tion 23.1 becomes now 


wo’ (x, h) = ¢’ (x,y, h) — F(a,h)). 
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Since 22.1 must hold for all values of c—y- a it must hold for a9 
and y = 0 so that we have 


¢’ (2, h) ¢’(h, 


and 


23.7 F(a, h)) = o(h, F (a,2)) 


and this must hold for every normal vector a. The formula 23.7 gives us 
the required algebraic integrability condition in the case now under consider- 
ation; relation 22.5 is, of course, a special case of this, F(a, h) becoming s(h) 
in the case when there exists a unique unit normal vector. 

Using the linear vector function f(z) instead of the tensor of the second 
rank ¢ we may write our relation as 


f(x) .F(h,a) =f(h). F(z, a) 
or, with the aid of 23.5, 

Since this must hold for an arbitrary a, 

s[f(x),h] =s[f(h), 


The generalization of 23.7 for the case of a tensor of higher rank is 
immediate. For rank three we have, for instance, outside of the differential 
condition the two algebraic conditions 


o[z, F(a, h),z] = F(a, x), 2] 
an 
23.81 y, F(a, h)] = o[h,y, F(a, 2) ]. 


§ 24. The Tensor Area. 


In higher spaces the Riemann tensor plays in some respects the same réle 
that the indicator plays in the case of an ordinary surface (and, more gen- 
erally, in the case N—-n=1). We considered the indicator as giving the 
infinitesimal change of the unit normal vector which corresponds to an infini- 
tesimal displacement along a curve on the surface, and its integral gave us 
the finite change of the unit normal vector which corresponds to a finite dis- 
placement along a curve. In the case which we are about to consider now 
the réle played by the curve will be played by a surface belonging to the n-S. 
considered ; instead of the change of the normal unit vector between the two 
ends of a finite curve (which constitute the frontier of that curve) we shall 
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consider the change in the orientation of normal spaces which corresponds 
to the frontier of a finite piece of surface, i. e. to a contour, and we shall see 
that in the case of an infinitesimal contour the corresponding infinitesimal 
change of the normal spaces is connected with the Riemann tensor. Before we 
come to this, however, we have to introduce a new notion which has the same 
relation to the frontier of a piece of surface as the vector chord, i. e. the vector 
joining the ends of an arc, has to these two points which constitute the fron- 
tier of this arc. We consider here an N.-E.S. 

The vector which joins the points A and B we agreed to represent by 
B— A; we know also that the consideration of a vector is equivalent to the 
consideration of its product with a variable vector—a linear scalar function 
(First paper, 3.4). In our case this tensor is (B—A).z, but it is more 
convenient to write it in the more complicated form 


1 (A—P).z 
weit 
where P is an arbitrary point on which the expression obviously does not 
depend. 
It is suggested, if we are given three points A, B, C, to consider the 

expression 

1 (A—P).2 (A—P).y 
24.1 B(A, B,C; =4 1 (B—P).« (B—P).y 

1 (C—P).x« (C—P).y 


This expression also does not depend on P; it is for given A, B, C a tensor 
of the second rank which is alternating or antisymmetric in the variables 
y, 


In the case of a plane, B(x,y) is simply the area of the triangle ABC 
multiplied by the so-called alternating product of the vectors x and y; in the 
case when NV = 83, i. e. in ordinary space, B(A, B,C;2z,y) can be presented 
in the form of the scalar product of a certain vector v by the vector product 
of x and y; this vector v is perpendicular to the triangle ABC and its length 
is equal to the area of this triangle; it is called the vector area of the triangle 
or a bivector; the vector area of a triangle as well as its generalization for 
any contour have been considered by Peano, Laisant and Fréchet (see the 
latter’s paper “Sur une généralisation des notions d’aire et de plan,” Nouv. 
An. (4), 4 (1904), pp. 241-248); the present writer used it in connection 
with integration in space (see § 27 below). 
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When WN > 3 it is impossible to get rid of the variable vectors x and y 
as in the cases N = 2 and N =83 so that we have to consider the thing as a 
function of z and y. Schouten uses the term bivector also in this case, 
(Compare e. g. J. A. Schouten and D. J. Struik, “ Einfiihrung in die neu- 
eren Methoden der Differentialgeometrie,” I, §1.) It should be noticed that 
making P = A we obtain a simpler but less symmetrical expression for 8, viz, 


(B—A).x (B—A).y 


24.11 B(A, B,C; =4 (C—A).2 (C—A).y 


If we are given a rectilinear contour A, B, C, D,--:, K, L or a number 


of points in a given cyclic order, we form 


24.2 B(A, B, C,---, K, L, A; 2, y) = B(A, B, P;2,y) 
+ B(B, C,P;2,y) + B(C, D, P; 
+ B(K,L,P; + B(L,4,P; 2,9), 


an expression which, obviously, does not depend on P; we call this the tensor 
area of the contour ABCD :--- KL. Given now a curvilinear contour I, we 
may approximate it by rectilinear contours and define the tensor area of I, 
which we shall denote by B(T; x,y), as the limit of the tensor areas of these 
approximation contours when this limit exists, which happens in a very gen- 
eral case. 

The notion may be extended by starting with four, or any number of 
points, instead of three; we have for the “tensor volume” of a tetrahedron 
ABCD the expression 

(A—P).« (A—P).y (A—P).z 
(B—P).« (B—P).y (B—P).2z 
(C—P).2 (C—P).y (C—P).z 
(D—P).« (D—P).y (D—P).z 

If we have any number of points (> 4) in a given spherical order (as 

analogous to a cyclic order) we may form the tensor volume for such a poly- 


hedron and by a limit process extend the notion to apply to a closed surface 


ete. 


§ 25. Contour variation of normals and a new introduction of the 
Riemann tensor. 


We introduced the indicator with the aid of the change b—a of the 
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normal unit vector; instead of the vector b —a we may here again write the 
tensor of the first rank (6 —a).z or 


| 1 @.2 
1 


25.1 


If N—n=r>1 this cannot be generalized immediately. We denote 
as before the tangent bundle at a point A by the same capital letter A, and 
the normal bundle we denote by the corresponding small letter a. The pro- 
jection of a vector z on the normal bundle we denote by affixing this letter a: 
thus %.; in the case r—1 this would be a(a.r). The expression 25.1 sug- 
gests the formation of 


1 La Ya 
25.2 5(A, B, C; y) 1 Yd 


when we are given three points A, B and C; this is a tensor of the N.-E.S. 
which characterizes the relative positions of the normal r-dimensional Eu- 
clidean spaces to one curved n-S. at the three points A, B, and C. It is, 
like the tensor area, an alternating tensor and, exactly as in the case of the 
tensor area, the definition of this tensor, which we call the variation of normals 
(normal Euclidean spaces or normal bundles) may be extended to apply to 
any rectilinear, and then to a curvilinear contour by first forming the 
expression 


B, P; x,y) + 8(B,C,P;2,y) +--+ 
and then passing to the limit. 


We thus obtain for a contour [ a tensor 8([;2,y) the variation of 
normals to the n-S. for the contour T. This corresponds to the finite change 
of the unit normal vector between two points. 

Now, according to our program, we pass to the infinitesimal variation, 
which corresponds to an infinitesimal contour. We begin with a triangular 
contour A, B, C considering A as a fixed point and B and C as variable points 


—A — ha, —k,s. We write 


approaching A in such a way that 
25.2 in the form 
Ly — La Yo — Ya 


If x and y are vectors of the tangent bundle at A we have %2_—0, ya =0; 


if 
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we may write in this case 74 for z and ya for y; 2» will be 4 — 24.2 ete, 
so that 
8(A, B, C; ya) 


TA.B— VA YA.B— YA 
TA.0 — VA 


Dividing both parts by Aw and passing to the limit we have (§ 17, Second 
paper) on the right hand side the Riemann tensor at the point A. A con- 
nection is thus established between the infinitesimal variation of normals and 
the Riemann tensor. 

But this is not the whole story; the Riemann tensor gives us only those 
values of the tensor expressing the infinitesimal variation which correspond 
to the values of x and y belonging to the curved space. In order to know 
the infinitesimal variation of the normals completely we have to consider not 
only 8(24, ya) but 8(2,y) for all values x,y. If we write in the general case 


2%, + Ya, We have 


Lp — La = — + Lad — Xa, 


where 2a» denotes the projection of the normal component (with respect to the 


point A) 2 of z on the normal bundle at B; but then 2, —%a = — (2u)B 
and we have (23.6) 


ik 
id 


lim 


i 


The full expression for the infinitesimal variation is therefore 
8(Za,ha) + F(ta,ha) (ta, ka) + F (aa, ke) | 
8(ya, ha) + F (Ya, ha) 8(ya, ka) + EF (ya, ks) 


Since s and § are perpendicular to each other, s being a normal vector 
and Fa tangential one (see 23.5), the products of the type s. F vanish and 


we only have 


s(ta,ha) 8(2a, ka) | F (taba) F(a, ka) | 
s(ya,ha) s(ya, ka) F(ya,ha) ka) 


The first term is the Riemann tensor. We see thus that the Riemann 
tensor gives us only a part of the infinitestimal variation of normals; we may 
say that it gives all the information concerning the variation which can be 
expressed in terms of the curved space, of its internal properties, because the 
other part involves vectors 2a, ya which require the consideration of the con- 
taining space, i. e. which refer to external properties (compare § 17, Second 


25.4 


paper). 
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Here we have a radical difference between the Riemann ‘tensor in general 
space and the indicator of a surface: the latter gives us full information 
concerning the change of the unit normal vector because this change can be 
controlled from the surface; the variation of normals cannot be controlled 
from the space (in general) and the Riemann tensor, which is an internal 
thing, only gives us partial information. 

But there are cases where the whole variation can be given by the Rie- 
mann tensor; this happens obviously, when the second part of 25.4 vanishes, 
i.e. for spaces for which we have 


25.5 F (Xa,h). F (Ya, k) = F (Xa, F(Ya, h) 


in these cases, we may say, the whole infinitesimal variation lies in the curved 
space and is, therefore, given by the Riemann tensor. Only in these cases can 

_ we hope to obtain the finite variation of normals, which corresponds to a finite 
contour, by integrating the Riemann tensor in analogy to what we had in the 
case of a surface. 

We shall not investigate here the peculiarities of spaces for which con- 
dition 25.5 is satisfied ; this would lead us to considerations which lie outside 
the scope of the present paper; in dealing with surface integrals we will. 
however, again be led to relation 25.5. 

We want now to get rid of the restriction which we imposed upon our- 
selves by considering only triangular infinitesimal contours. Let us consider 
a rectilinear infinitesimal contour, constituted by the points A +h, A+k, 
A+lA+m,:::,A+n,A-+h. It is clear from the preceding that the 
corresponding variation of normals is given by p(z,y;h,k) + p(z,y; k,l) 
+ p(z,y31,m) ++ +p(ay;n,h). 

From the expression 17.2 for the Riemann tensor we conclude that, for 
fixed z and y, p(z,y;h,k) is a sum of terms of the type 


ph pk 
gh qk 


where p and q are vectors belonging to the bundle A. The infinitesimal vari- 
ation of normals which we are considering has, therefore, for its expression 


{ 


where the summation is extended over a certain number of pairs of vector 
Pr, G1» Jz ete. 


pn ph 
qn gh 


ph pk 
gh qk 


qk ql ql qm 
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Now the tensor area of the contour A+h, A+hk,---, A+nisg 


«xk zk al al zn sh 
yh [ye yt | yt ym] yn yh 


25.7 | 


Suppose that we take another contour with the same tensor area; this 
means that 25.7 has for the new values of h, k, 1,- - -, m, n and for arbitrary 
values of 2, y, the same value as for the original h, k,1,---, m,n. Making 
x, y in turn equal to p,, 91; Po, J2 etc. we see that all the terms of 25.6 do 
not change when we substitute new values of h, k, l,- - +, m, n for the old 
ones, provided the new contour has the same tensor area. The infinitesimal 
change of normals for an infinitesimal contour T in the neighborhood of a 
point A of the n-S. 1s thus seen to depend solely on the tensor area of this 
contour; and of course it is a linear function of this tensor area. This com- 
pletes the analogy with the indicator considered as a linear function of the 
displacement vector giving the infinitesimal change of the unit normal vector, 

We may add that this gives a plausible interpretation of the Riemann 
tensor as a bivector tensor as Schouten puts it (loco citato, p. 34). 


§ 26. Contraction. 


In order to deal with integrals extended over surfaces, volumes etc. we 
must now introduce the operation of contraction, the use of which we have 
been able to avoid hitherto. 

We start with the identity 


26.1 ty (tp2) 
v=1 

where 1;, t2,° * *, t, are n mutually perpendicular unit vectors of an n-dimen- 
sional bundle and which holds for every vector x of the bundle. We shall 
have occasion to use expressions of this type, i. e. sums of n terms each term 
of which contains twice one of the n mutually perpendicular unit vectors; 
if we have a tensor or a multilinear vector function $(2, y,z,° - -), the ex- 
pression which we have in mind is 


n 
=X $ (ty, 
v=1 


It is easy to show that such an expression does not depend upon the particular 
set of n mutually perpendicular unit vectors which we use, i. e. that it does 
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not change if we substitute one such set for another. We therefore introduce 
for such a sum the notation ¢(J, J, z --~-) replacing each of the two vari- 
able vectors for which 1, stand in the sum by a Roman numeral. If such an 
operation has to be performed more than once on the same tensor we use each 
time another numeral, so that #(J, IJ, I, IJ) stands for 


This operation is essentially the same as the operation of contraction 
(Verjiingung) of tensors in the usual form of tensor calculus with the dif- 
ference that we have here no distinction between lower and upper indices. 
The Roman numerals stand here for the dummy suffices. 

With the notation introduced, identity 26.1 takes the form 


26.11 I (Iz) 
' multiplying this by y we obtain 
26.2 ry = (Ix) 


which means that the product of two vectors is the sum of the products of 


their components. 
If we have a scalar linear function ¢(x)—a tensor of rank one—we may 


write 


26.3 = (Iz) ] = Uz) = [¢(2) 


The expression ¢(J)I is a vector which written explicitly is 


i, 


if we denote this vector by g, we have ¢(z) —g.z. The operation of con- 
traction permits us thus to write explicitly in terms of a scalar linear function 


the vector which is the constant factor of it (compare property 1 § 3, First 
paper). 

It may be said that our definition of contraction contains traces of a 
coordinate system; this may be avoided to a certain extent by using the 
vector product of n—1 vectors but the expression, though more intrinsic, 
becomes more cumbersome and it was not thought worth while to introduce 
it here. We only shall mention that in the case of 3 dimensions we have 


$(I, 1). (xyz) = $(2,4°2) +4(y,2°2) 


where a*y denotes the vector product and zyz the triple scalar product, and 
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x, y, 2 are any three not co-planar vectors and in the case of four dimensions 
we have similarly, 


$(I, 1). (xyzu) = ¢(xyz, vu) —o(yzu, x) + y) 2) 


where xyz is the vector product and zyzu the quadruple scalar product. 
The operation of contraction may be used to pass from one form of the 
indicatrix to another; for instance (23.5) 


F(a, 2) =y(2,I,a).I = [a.s(2,1)].1. 


We are not going to enter into the details concerning the products men- 
tioned above; the vector algebra in 4 as well as in 3 dimensions can be devel- 
oped intrinsically (see e.g. a paper by A. Mochoolsky in the Memoirs of the 
Research Institute, Odessa, February, 1924 (in Russian) ) ; here for the sake 
of completeness we only shall state that, in three dimensions the vector product 
2° y is the vector whose components are the determinants of the matrix 


Le 


Yo Y3 


the triple scalar product xyz in three dimensions is the three row determinant 


formed with the components of x, y, z; the vector product xyz in four dimen- 
sions is the vector whose components are determinants of the matrix 


Zo 23 24 


and the scalar quadruple product is the four-row determinant of the com- 
ponents of four vectors. 


$27. Surface Integrals in Curved Space. 


If we have in a n-S. an alternating (antisymmetric) tensor field of the 
second rank ¢(A,72a,ya) and a surface 3, we may integrate this field over 
the surface. If we imagine the surface divided into “ infinitesimal triangles” 
so that two of the sides of each triangle may be considered as displacements 
of one of its vertices A, we form the sum of expressions ¢(A,dA,5A) and 
this sum might be considered as the surface integral. But it is more satis- 
factory to consider any finite division of the surface by means of a net of 
curves on it. We form the tensor area B(T,; z,y) for the contour which 
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ions constitutes the frontier of each part, we multiply this by the value $(Ay, 2, y) 
of the given tensor field at some point A, belonging to the same part of the 
surface and contract this product so as to obtain 


$(Ay, LIT). B(T,; 1,11); 


the we form the sum of these expressions for all the parts into which the surface 
was divided and consider the integral as the limit (if any) to which these 
sums tend when the diameters of the parts (i. e. the upper limits of the 
distances between two points belonging to the same part) tend toward zero. 


n- It is not the purpose of this paper to enter into details concerning the 
el- definition of integrals and the conditions under which they exist; we are in- j 
he terested here in questions of a less subtle nature; it might be interesting, how- : 
ke ever, to notice that the above definition of a surface integral does not suppose 
ict the existence of tangent planes to the surface. Integration in ordinary space 


was considered from a similar point of view in a paper by the present writer 
which appeared (in Russian) in the Bulletin of the Physico-Mathematical 
Society of Kazan in 1915. 

If we consider the surface in a Huclidean space and have in this space 


nt a tensor field of a rank higher than two, which is alternating in two argu- 
ments, say (A, 2, y,2,° °°), with 
$(A,2,y,%,° =— *), 


we may integrate it over the surface by giving to the “superfluous” argu- 
ments z etc. fixed values. But if we consider the surface in a curved n-dimen- 
sional space this is not possible and we have to resort to the device, used in 
- §21, of giving to z and other arguments values which are projections of 
fixed vectors in the containing NV.-E.S. We obtain thus as the result of in- 
tegration of a tensor field of rank k (which is alternating in two of its argu- 
ments) over a surface a tensor of the containing space of rank k — 2. 
Integration over a volume, i.e. a spread of 3 dimensions presents after 
this no difficulties: the integrand must be a tensor field alternating in three 
of its arguments; we divide the volume in parts; combine the tensor volume 
: of the frontier of each part with the value of the integrand at some point 
) of this part, by multiplying and then contracting three times; we give to 
superfluous arguments values which are projections of fixed vectors of the 
containing space; form the sum of the numbers so obtained; and pass to the i 
limit making the diameters of the parts tend toward zero. The integral is a | 
tensor of rank & — 3 of the containing space. i 
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§ 28. Integrability Conditions for Surface Integrals. 


We shall take over from the usual coordinate exposition, similarly to what 
was done in § 22, the integrability condition for the surface integral of , 
tensor field of rank two. This condition, viz., the condition that an integn| 
of ¢(A,24,ya) over a surface depends only on the contour of the surface, 
i. e. that it is the same for two surfaces having for their frontiers, the same 
contour provided that one can be continuously transformed into the othe 
without leaving the part of the n-S. in which $(A, 2a, ya) is. considered, is 
that in this part ¢ satisfies the relation 


28.1 $’(x,y,h) + $’(y,h, x) + $’(h, =0, 
where y, 2) is the differential of ¢(z, 


If we are given a tensor field of higher rank than two, let us say, for 
example, a tensor field of rank four ¢(z,y,z,u) which is alternating in the 
arguments z, y, and we integrate it with respect to these arguments using 
the method indicated in § 27 we have to write that the integrability condi- 
tions are satisfied for any tensor field w(z, y) which is obtained from 
$(z,y,2,u) with the aid of two vectors c,q of the N.-E.S. which contains 
the n-S. in which the tensor field ¢ and the surfaces over which we integrate 
are given, Viz., 


o(A, TA; Ya) (A, YA, CA; ga). 


Just as in § 23 we obtain the conditions 
F(a,h),u] + oly, h, F(a,2),u] + 2, F(a,y),u] =0 
an 
28.3 ¢[z,y,2,F(a,h)] + oly, h, 2, + 2,2, F(a,y)] =0 


which must be satisfied for all values of the normal vector a in addition to 
the differential condition 


28.4 $’(z,y,2,u,h) + $’(y,h,2,u,x) + =0 
where ¢’(z,y,2,u,h) is the differential of y, z,u). 


Let us apply this to the Riemann tensor p(z,y;z,u). The differential 
condition is always satisfied; in fact it is nothing but the Bianchi relation 
12.5 (First paper) proved to hold in the general case toward the end of § 18 
(Second paper). It is worth while to remember here that the differential 
condition for the indicator of a surface coincides with the Codazzi relation 
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which thus plays for the indicator the same réle as the Bianchi relation plays 
for the Riemann tensor. We need only consider one of the algebraic condi- 
tions 28.2 and 28.3 since the variables z and wu enter in p in the same way 
(only the sign is changed when they are interchanged so that 28.2 is equiva- 
lent to 28.3 in this case). If we use expression 17.3 (Second paper) for 
the Riemann tensor, 28.2 becomes 


s[y, F(a,x)] s(y,u) 
s[h, F(a,z)] s(h,u) 


s[h, F(a,y)] s(h;u) 
F(4,y)] s(z,u) 


s[2, F(a,h)] s(z,u) 


F(a,h)] s(y,u) 


Expanding the determinants we obtain a sum of six products, the first of 
which is s[z, F(a, h)].s(y,w); we may identify this with the last term of 
23.5 considering s(y,u) as a, x as ka, and Y(a,h) as ha; using then the 
form of the middle term of identity 23.5 we may write our product as 
F(a,h). F¥[s(y,u), x]; applying this to all the six products of the expansion 
of 28.4 and changing the order of terms we may write it as 


F(a,h). F¥[s(y,u), 27] — F(a, 2). F[s(y,u),h] + F(a, zr). F[s(h, u), y] 
— F(a,y).Fls(h,u),c] + F(a,y).F u), h] 
— (a,h).F[s(z,u),y] =0. 


This vanishes when condition 25.5 is satisfied; this is very natural because in 
this case the Riemann tensor gives the full infinitesimal variation of normals 
and therefore its integral over a piece of surface gives the finite variation of 
normals for the frontier of this piece, i.e. for the contour which limits this 
piece and, of course, this does not depend upon the surface—just as the 
integral of the indicator along a curve depends only on its end points. In 
the general case the Riemann tensor does not give the full infinitesimal varia- 
tion, its integral does not give therefore the finite variation and there is no 
reason why it should be determined by the contour alone. 

As our next example we take the contracted Riemann tensor which we 
write, according to § 26, as p(a,I;z,I). This is a symmetric tensor of rank 
two which we denote by e(z,z); but instead of considering it as such, we 
make of it a tensor of rank n: we substitute for one of its arguments, say z, 
the vector product of nm —1 vectors. If we take, for instance, n = 4 we have 
to make z= uvw, i. e. the vector product of the vectors u, v, w (see end 
of § 26). 

Now we have in e(z, uvw) a tensor which is alternating in the 3 argu- 
ments u,v, w. It may be therefore integrated over a volume (3-dimensional) 
and gives a vector of the containing space; we shall see what are the integra- 
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bility conditions, i. e. under what conditions the integral vanishes over every 
closed volume. If we write 


y(u,v,w) = (C4, U,V, W) 
the conditions in question are given by 
y/ (u,v, w,h) —y(v,w,h,u) (u,h, u,v) —y(h, u,v, vw) =0. 
This again gives a differential condition for e viz., 
(x, uvw, h) —e(z, vwh, u) + whu, v) huv, w) =0 
which may be written as 


1,1) =0 


and an algebraic condition, viz., 


«| F(a,h), —e [ F(a, v), vwh] + ¢ (F(a, v), whu] —e [F(a, w), huv] =0 


or 


e[F(a,I), I] =0. 


The discussion of spaces where these equations are satisfied will be given 
_ elsewhere. 


§ 29. Conclusion. Integral Expressions Which do Not Require the 
Consideration of the Containing Space. 


We had to go out of the space in order to consider integrals in the 
preceding sections as soon as the integrand had superfluous arguments, i.e. 
more than one in the case of a line integral, more than two (in which the 
tensor is antisymmetric) in the case of a surface integral, and so on. But 
it is desirable to be able to form in such cases also integral expressions not 
involving the conideration of the containing space. In some cases it is pos- 
sible. Let us begin again with line integrals and let us take as the integrand 


a tensor field of the second rank which is symmetric in its arguments, so that 


we have 
29.1 $(z,y) = o(y, 2). 


It can be proved that such a tensor can be presented in the form 


29.2 = + (ax) (ay) + (br) (by) (cx) (cy) +° 


where a, b, c,+ + + are mutually perpendicular vectors, the number of terms 
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being equal to the number of dimensions of the space and the + signs being 
used in order to avoid the introduction of imaginary vectors. Now the pro- 
ducts az, bx, cx, etc. are tensors of the first rank; they can be integrated, 
therefore, without leaving the curved space; we introduce the notations 


29.3 P= fa.dA, Q= f b.dA, R= f c.dA, 
and form the expression 
29.4 


using the same signs as in the corresponding places in 29.2. The expression 
29.4 is the sought integral expression. Let us consider as an example the 
case of an ordinary (two-dimensional) surface and as the integrand let us 
take the indicator; if 1, 7 denote the two unit vectors which belong to the 
principal directions and A, » are the curvatures we have a= iVa, b =j Vp 
so that it will be convenient to introduce the notations VA=a, Vu=—=B. 
Let us ask ourselves when the integrals P and Q are independent of the 
path. Since the integrand of P is w(x) —a.(ix) the condition 22.1 takes 
the form 


29.51 a’(h). (tr) + a[i’(h) = a’(x). (th) (2) .h]; 
for Q we have 


29.52  B’(h). (7x) + = B’(x). (jh) + 


On the other hand the Codazzi equation for the indicator, which now has” 


the form 
a’ (1x) (ty) + (jx) (jy), 
gives 
2a.a’(h) (1x) (ty) + a*[0’(h) (ty) + (iw) (h) -y] 
29.6 + 28.B’(h) (jx) (Jy) + (h) + (92) (2) 
2a.a’(x). (th) (ty) + (x) .h] (ty) + (th) -y] 
+ 28.8’ (x). (jh) (jy) + -h] (iy) + (9h) [7’ (2) -y]. 


The simultaneous existence of equations 29.51, 29.52 and 29.6 leads to 
the relation (see for the method used the author’s paper in the January 
number of the Transactions of the American Math. Society for 1925, p. 123), 


a’? + or 


so that only for a minimal surface are the integrals P and Q, and therefore 
the expression 29.4, independent of the path. These questions connected 
with minimal surfaces will be treated in detail elsewhere. 
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Let us consider now the case of a surface integral and a tensor field of 
rank four ¢(2,y; u,v) with the property 


u,v) = (u,v; 


as the integrand in a 4-S. Suppose that @ may be presented in the form 


ku ky ku kv 
lr ly lu Ww 


is wm 4 
je ju 
Precisely such a tensor occurs in physics, viz. the second part of the Riemann 
tensor of space-time. It is natural, in order to find an integral expression 


to put again VA—a, Vi= B and to consider the expression 


P+ 


where 


P= fa 


i.dA i.8A 
j.dA 


k.dA &k.8A 


This expression has in certain cases a physical meaning, being connected 
with the square of the electric charge; in other cases it is possible to obtain 
by a slight modification another more general integral expression of the same 
type. We shall not go into details concerning this type of integral expres- 
sions. The purpose of this section was only to indicate by way of conclusion, 
that the integrals involving the consideration of vectors of the containing 
space are not the only integral expressions which can be formed with tensor 
fields having superfluous arguments. 


BALTIMORE, JUNE 14, 1925. 
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On Certain Symmetric Sums of Determinants. 
By L. DInzs. 


The first theorem of this paper has to do with a certain transformation 


of determinants which may be defined as follows. If D is a determinant of 


order n, and if a,, a2, ***, an is any sequence of numbers (not necessarily 
ann distinct) each chosen from the set 0, 1, 2,---, n, then by the notation 
we shall mean the determinant of nth order obtained from D by shifting / 
each element of the first column up a, places (thus eliminating the upper a; : 
of them from the determinant), shifting each element of the second column qi 
up a, places, and so on to the nth column, and then filling in the vacant 1 
places at the bottom by zeros. From the definition it is obvious that 4 
‘a To ..oD =D, while if every a; is positive or if any a; is equal to n, then i 
Daa,..an D = 0. Between these extremes are the cases of interest. 
a For a given sequence a, a2, ***, an, we shall consider all the distinguish- 
ee able * permutations of the sequence. 
By the notation 
| 
or we shall mean the sum of all the determinants Taj,a;, -- ai,D where ai,, ai,,°*, 


a;, is a permutation of a;, a:,°*, an. The theorem? to be proved in §1 is 
then to the effect that if at least one element of the sequence ai, a2,**, an 18 
positive 

ST . an D = 0. 


For example for n = 3, we have 


As bz a; be as by as 
ST 001 b; be Ds = bo C3 Co Ds Ci be Ds = 0, 
Ci Ce Cs Ci Ce 0 Cc, Cs 0 Co Cs 


*If no two elements of the sequence are equal then all its permutations are 
distinguishable. If certain elements of the sequence are equal, permutations which 
differ only by the interchange of equal elements are not distinguishable. In this case 
only one representative of the set of non-distinguishable permutations will be con- 
sidered. 

¢ Theorems somewhat similar to this one and to its generalization stated in §1 
have been given by Sylvester, Hammond, and Deruyts. See Muir’s History of Deter- 
minants, Vol. IV, pages 2 and 25. ; 
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and other vanishing sums corresponding to the sequences (0, 0, 2) ; (0, 1, 1); 
and (0, 1, 2). 

The theorem of § 2 has to do with sums of determinants of nth order 
from a persymmetric * array of n rows, the elements however being affected 
by binomial coefficients * according to a definite law. 
For n = 3 the theorem would have reference to a persymmetric array 


Denoting the binomial coefficients by the notation 


1-2: -k 


(m)k 


we may state a simple case of the theorem as follows: 


(J-1) (k+2) (7+-3) 
= (J7+2) 2C5+2 (k+3) (1+-4) 


 Clag 


= 


jtk+l=p 


Class 


the identity holding for every non-negative integer p. 


Theorem: If D is a determinant of order n, and a, a2, ***, an 18 @ 
sequence of numbers (not necessarily all distinct) each chosen from the set 
0,1,---, n, and at least one of them positive, then 


BT eas. Gn Dz 0O. 


For n = 2 the theorem is obviously true since 


G11 G11 


ST 
01 


0 


Assuming the theorem to be true for n —r—1, we will prove it true 


*For literature on persymmetric determinants and on determinants whose ele- 
ments are binomial coefficients, see Muir’s History of Determinants, Vol. IV, Chapter 
XIII, and Vol. III, Chapter XX. 


it 
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ay 


for n= rT, thus establishing its general validity by mathematical induction. 
If we denote by D‘‘) the algebraic complement of the element ai; in 
the determinant 


we find upon consideration of the expansion of D according to elements of 
its first row, that 


= r 


k=1 j=1 


the notation >’ indicating that in case a; has the same value for two or 
k=1 


more values of k, the summation as to & shall include only one term repre- 
senting such values for each value of j. 

Now with respect to the symmetric sum S on the right side of (1), 
there are two possibilities. First, if any one of the numbers aj, a2, °*, axa, 
ans, °*, ar is different from zero, then the symmetric sum vanishes on account 
of our assumption that our theorem is true for determinants of order r—1. 
On the other hand, if a; = = = = =a, = 0 then by hy- 
pothesis a, 540. And for the value of & which gives rise to this situation, 
the right side of (1) reduces to 


j=1 


which vanishes by an elementary theorem on determinants. Hence the entire 
right side of (1) vanishes, and our theorem is established. 

A number of variations of the theorem immediately suggest themselves, 
For example, the “shifting” involved in the transformation 7 might obvi- 
ously be made downward, to the right, or to the left, as well as upward. We 
will state one extension which is not quite so obvious, and which will be 
used in § 2. 

With each element a, of the sequence a;, as, -**, an we will associate a 
sequence of n — ax elements 


subject only to the restriction that if a; —a,;, then s;==s;. And we will 
denote by 
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the determinant obtained from D by shifting the elements of the kth column 
up a, places, filling in the vacant places by zeros, and then multiplying the 
upper n —ax elements of the resulting column by the respective elements of 
' the sequence s;, this procedure to be followed for each column of D. 


By the notation 
ST 4,83 8n D 


An 


we shall mean the sum of all determinants 


* 


where ai,, ai,, °°, ai, is a permutation of a;, a2, °*, an. We may then state 
the generalization in the form: 


If D is a determinant of order n, and ay, a2, **, an is & sequence of 
numbers chosen from the set 0, 1, 2,-°--, , of which not all elements are 
zero, ANd $1, So, **, Sn are sequences of form (2); then 


D=0. 


The proof may be made by an induction similar to that used in the simpler 
theorem. 


§ 2. 

In this section, we suppose the elements of the persymmetric array 
C1 Co Cg 
C2 Cy C4 


Cn Cn+1 Cns2 


to be well defined for as many columns as the statement of the theorem 


requires. 
Denoting the binomial coefficients by the notation 


m(m—1)---(m—k-+1). 
(m)c= 


we will prove the following 


jit 


| 
| 


mn 
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Theorem: For any non-negative integer p, and any sequence of non- 
negative integers 11, 


) (jo+re+1 ) (jntra+1 ) jntratl 


For the proof we shall make use of the following easily deducible relation 
between binomial coefficients : 


(j-+r+h) ran 


= (j+r+n) ren—(n—h) 1 2—1) 2 ran 
— 


which holds for h —1, 2,---, n. 
By means of this relation, the elements of any column in (3) may be 
replaced respectively by polynomials 


= (jATrt+1) 

(j+17+2) ranCjarse 


(j-+r-+n) rinCj+r+n-1— ( 1 ) 1 (j+r+n—1 ) renCj+r+n-1 
(j-+r-+n) 
the symbols j and r carrying subscripts indicative of the particular column 


considered. 
Each determinant of the sum (3) can therefore be replaced by the sum 
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of n” constituent determinants (many of which obviously vanish) with the 
property that in each column all the elements contain a common binomial co- 
efficient factor. The typical form of these columns will be 


(—1)*(m—1) 
(—1)*(n—2) 
(4) 


(—1 ) (a) a (j-+r-+n—a) ranCj+rsn-a 
0 


0 

the last a elements being zero. There will be such columns corresponding to 
a= 0,1,--:,2—1, and we will say that the column (4) has the character a. 

Consider now any set of numbers a, as, ***, an, each chosen from the 
sequence 0, 1, ---, n—1, not necessarily distinct, but at least one of them 
being zero. Consider the constituent determinant in which the first column 
(involving the index 7,) has the character a;, the second column (involving 
the index j.) has the character az, and so on. We will say that this con- 
stituent determinant has the character (a;, a2, ***, an). Out of this deter- 
minant we can take the factor 


(jn Tn -++n—an) 


leaving a residual determinant of which a typical column is 


(n—1 ) aCj+r+1 
(n—2 ) aCj+r+e 


a Cisrsn-a 


0 


0 


If (ai,, ***, is any permutation of (a;, a2, an) there will 


Dines: On Certain Symmetric Sums of Determinants. 255 


be a constituent determinant having the character (ai,, ai,,°**, a:,), and out 
of it we can take the factor 


(jn+rn-+n—ai, ) 


which in general will be different from (5). 

It now becomes necessary for us to note a significant property of the 
sum (3) which we have not previously stressed, and which will enable us to 
correlate the products (5) and (6). The notation 


. . +jn=p 
indicates that the summation is to be made for every combination of indices 
jy joy °° * > Jn Of which the sum is p. Hence if in a determinant affected by 
this summation sign we make a transformation of indices 


je = J'2 + jn = J'n + Qn, 


there will be no change in the content of the summation, provided a, + a. + 
an = 0. 
If therefore in the determinant with character (ai,, ai,,-°*, ai,) we make 


the substitution 


the content of the summation will not be affected. And upon dropping the 
primes from the indices 7’;, j’2,° * +, j’n in the product (6) after transforma- 
tion it is seen to be identical with the product (5). 

By means of similar transformations we are able to reduce each deter- 
minant whose character is a permutation of (a, a2,° °°, an) to a-form such 
that the identical factor (5) can be taken from it. 

It is not very difficult to see further that the residual determinants are 


obtainable by a transformation 7'.,s,..;, from one and the same determinant. 


We may in fact write the sum (as indicated by = ) of all the con- 


+jn=p 
stituent determinants whose characters are permutations of (a;, a2, *, an) 


in the form 


(7) a (—1) +40 (9, +-n—a1) ran rein’ * 


. +in=p 
ST — D, 
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4 C jntrntn-an 
and 

Se = ((m—1) a,» *, (ax)ax)- 
We are now ready to draw our conclusion. If the set (a:, a2, --*, ag) 
contains an element different from zero, the sum (7) vanishes by the general. 


ized theorem of §1. Therefore the only non-vanishing part of the sum (8) 
is given by (7) for a, =a,—=-*+-—an,=—0. This value is evidently 


+jn=p 


Cigergen Cintrntn | 
and this is the value as stated in the theorem. 
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A Direct Solution of Systems of Linear Differ- 
ential Equations having Constant 
Coefficients.* 

By J. A. NyswANDER. 


1. Introduction. The theory of linear differential equations having 
constant coefficients began with Euler, about the middle of the eighteenth 
century, and was continued during the latter part of the century by other 
mathematicians, the most noteworthy of whom were Gauss, Laplace, and 
Lagrange. All of these writers except Lagrange treated a single equation of 
low order; Lagrange, however, was able to solve those special systems of 
linear differential equations in a few dependent variables for which the roots 
of the characteristic equation D(A) —0 were distinct. During the interval 
of about seventy years between Lagrange and Weierstrass, various symbolic 
methods were devised for solving a single equation of the nth order, while 
for the general simultaneous system of equations a treatment was given only 
for the simple case in which D(A) —0 has distinct roots. In this period 
special systems of equations for which D(A) —0 has multiple roots were 


solved by regarding them as the limiting cases of corresponding systems of 
equations for which D(A) —0 has distinct roots. This method is unsatis- 
factory logically and practically impossible of application. 

The first complete treatment of a system of differential equations in n 
dependent variables was given by Weierstrass.t He started with the system 


(1) = +° 4 } n), 
where 2’; == dz;/dt, and made linear transformations on certain bilinear forms 
obtained from (1), thus securing a simplified equivalent system of equations 


of the form 
1615 


*Read before the American Mathematical Society at Chicago, April 18, 1924. 


+ Weierstrass, Collected Works, Vol. II, pp. 74-75. 
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where A, is a root of order j of the equation D(A) —0. The properties of 
D(A) pertinent to the problem were expressed in terms of the elementary 
divisors associated with D(X) and its minors. Although Weierstrass’ process 
is logically sound it is indirect and tedious in application. Solutions of the 
problem were also given by Weber* and Frobenius.t Their methods are 
essentially the same in principle as that of Weierstrass and differ from it only 
in certain minor details. 

A method of treatment { which is entirely different in character from 
that of Weierstrass was given by Thomé. By differentiating a single one of 
the differential equations n —1 times and replacing the x’; at each step by 


S ax, Thomé arrived at a single differential equation of the nth order in 
coe of the z;, say z,. This procedure implies that a certain determinant of 
the coefficients of v2, --~- , % in the n differentiated equations is different 
from zero. Each of the n solutions of this differential equation of the nth 
order gives rise to a set of simultaneous linear equations in 22,- - - , &» which 
can be solved. If, however, this determinant is zero the procedure just indi- 
cated leads to a system of non-homogeneous differential equations which must 
be solved in the process of arriving at the solutions of system (1). Thomé’s 
treatment, although simpler in theory and application than any of the earlier 
methods, is likewise indirect and laborious to apply. 

The present paper treats the general system of differential equations by 
a method which is both elementary and direct, and which at the same time is 
much simpler of application than any of the previous methods. A general 
form for the solution of system (1) is obtained, from which, in any particular 
case, the desired solutions can be written at once in terms of the minors of 
the determinant D(d;) and of their A-derivatives. 


2. The Solutions when D(A) =0 has Distinct Roots. Before solving 
the problem in its greatest degree of generality the two cases for which the 
characteristic equation D(A) —0 has distinct roots, and a triple root, respect- 
ively, will be treated. 


Consider the simultaneous system of equations 


(1) x’; = + +. Aintn 


* H. Weber, Lehrbuch der Algebra, Vol. II, pp. 41-42. 
¢ Crelle’s Journal, Vol. 86, pp. 42-72; 146-209. 
¢ Crelle’s Journal, Vol. 131, pp. 8-24; Vol. 133, pp. 1-17. 
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where the aij are given constants. Let the question be raised whether or not 


they have a solution of the form 


(2) 
in which the A; and A are constants. On substituting (2) in (1), it is found 
that if (2) is a solution, then the A; and A must satisfy the conditions 


(ii + +++ + dinAn = 0, 


(3) 


Equations (3) are linear in the A; and will have a common solution in which 
not all A; are zero if, and only if, \ is a root of the characteristic equation 


Suppose that D(A) =0 has n distinct roots A1,° °°, An. For each of | 

these n roots equations (3) have a solution for which not all the A; are zero. 

Since for each 4 = A; at least one first minor of D(A;) does not vanish, one 

can express the n solutions of (3) in the form 

(5) (7,+=1, n), 

where the c; are arbitrary constants and the 8;,; are the first-minor cofactors 

in D(A;) resulting from the solution of (3) for the A; by determinants. In 

general, the second subscript & of 8jxi is used 4a indicate that 3jx4 Tepresents 

a kth minor cofactor (of order n —k) of the determinant D(;). 
On substituting (5) in (2), one obtains the n solutionsof system (1) } 


(i=—1,--- 


Aj = 


(6) 


Ln = . 

A proof of the fact that (6) form a fundamental set of solutions of (1) 
will be omitted here since a demonstration of this property for the general 
case will be given in Section 7. : | 

3.. Solutions of System (1) when D(A) =0 has a Triple Root. Let 
An be the n roots of D(A) = 0, A=A, being a triple root. The 
following cases, classified with respect to.the properties.of .D(A,), can arise: 


i 

| 
| 
| 
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All first and second-minor cofactors (8j:, and of D(A1) ar 

zero, while at least one third-minor cofactor (8;,x) is not zero, 
II. At least one 8;,, is not zero. 

III. All 8; are zero, at least one 8; is not zero, and at least one 8i1x(A) 

has a A-derivative that does not vanish for ’ = dj. 


Since the A-derivative of every (j—1)st minor of D(A) is a sum of 
jth minors of D(X), it is evident that 


d*D (A) = bis. 
Thus if A=), is a triple root of D(A) = 0, the case for which all 8:3, all 
Siox = all 313, = 0 cannot exist. It would appear that under III there would 
arise the sub-case obtained by replacing the last hypothesis of III by the one 

that every 


d 
E =0, 


[ ~0. 


That this cannot occur when A=), is only a triple root of D(A) =—0 is 
readily seen by considering D(A) when reduced to its normal form.* 


at least one 


CasE I. In Section 2 it was seen that if the system of differential 
equations 


has a solution of the form 
(2) = Aje™t (t=1,---,n), 


then the A; must satisfy the equations . 
n 

(3) (aii —A1) Ai + 5&4). 
k=1 


Since by the hypothesis of this case D(A,) is of rank n — 3, the solution of 
(3) may be written in the form 


(4) Ai = 8is1€1 + + 0), 


where ¢,, Cz and cs are arbitrary constants, and the 83, are the third-minor 
cofactors in D(A,) obtained in the solution of system (3) by determinants. 


* Bocher, Introduction to Higher Algebra, p. 267. 
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One observes that in consequence of the nature of solution (4), at least two 
of the (1 are identically zero for k = 1, 2, 3 respectively, 
while at least one of the 8;3x, in each of these three sets, is different from 
zero. On substituting (4) in (2), one obtains as the three solutions of equa- 


tions (1) 


(5) =m = (fam ++, n; 2, 8). 
The remaining n — 3 solutions of system (1) corresponding to 
An are given by (6), of Section 2, fori1—4,- - -, m respectively. 


Case II. At least one 81, is not zero. It will be shown that in this 
case the solutions of equations (1), associated with A = ),, have the form 


where the A,‘ are constants. On substituting (6) in (1) and equating 
coefficients of #7, ¢ and 2°, one obtains as sufficient conditions that (6) shall 
be a solution of (1) 


(7) (as — Ar) + = (i=1,---,n; ki), 
k=1 

(8) (aig — AL) + k At), 
k=1 


(9) (a4 —A1) As + = Ay ki), 
k=1 


It will now be shown that these conditions can be satisfied. 
In order to obtain the solutions of (7), (8), and (9), start with the 
identities 


(10) (A) + + An-1,n0n1 (A) = 0, 


where 8::(A),* *, 8n1(A) are the cofactors of the elements @n1, °° *, @nn—A 
respectively in D(A). The third subscript of the 8j:x, which indicates the 
value of A under consideration, will be omitted except when it is necessary to 
distinguish notationally between several groups of cofactors of the same order. 
In writing the identities (10) it has been assumed, for convenience, that the 
first nm —1 rows of D(A,) contain a non-vanishing first-minor. 
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On differentiating (10) twice with respect to A and setting A= Aj, ity 
found that the relations 


(11) — Bax + 2 = 0 


(12) (ais — Ai) + = Sir peas 
=1 


(13) (ait —A1) 841 + = L 


are satisfied. In relations (12) and (13) the notation 


ix ( ) | ik—= | Six | )\ 


has been used. Relations (11), (12), and (13) give, at once, the following 
solutions of equations (7), (8), and (9) respectively : 

(14) Ag™ == n), 

(15) A; + Codi1 1, n), 

(16) == + 62841 + 

where ¢,, C2, and cs are arbitrary constants. On substituting in (13) the 
values of the A;‘/ just determined, one obtains three solutions of (1) which § 
can be written in the following form: 


The remaining n— 3 solutions of (1) are given by (6), of Section 2. 


CasE III. All 8:1, are zero, at least one 8:2, is not zero, and at least one 
8’i1x is different from zero. 

By the procedure employed in Case I, one obtains the two purely expo- 
nential solutions of (1) 
(18) = + (t= 2). 


The remaining solution of System (1) corresponding to =A, is obtained 
from the first of the three solutions included in (17). For the hypotheses of 
Case III this solution becomes 


(19) = (8G + 


The remaining n — 3 solutions of (1) are given by (6) of Section 2. 
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4, Fundamental Sets of Solutions. For each of the cases considered 
in the two preceding articles n solutions of equations (1) were obtained, but 
it was not shown that they form a fundamental set. Although this property 
will be established in Section 7 for every case that can arise, yet it seems 
advisable at this point to illustrate the method by proving the independence 
of the n solutions in Case II, of Section 3. 

The n solutions of the system of differential equations (1) will form a 
fundamental set if, and only if, their determinant A(¢) is different from zero. 
Suppose that the solutions 


(20) 


do not form a fundamental set, and hence that A(t) —0. One easily veri- 


fies that 
where 
(21) A(0)= 


in consequence of the assumption A(t) = 0. 
The assumption (21) implies that there exist n constants ¢,°°,- - -, ¢n™® 
not all of which are zero, such that 


(Vj=—1, n). 


For this special set of values of the c;, equations (20) become 


(23) 
i=4 


Thus (23) is a solution of (1) having, from (22), the initial conditions 
2,(0) = 2,(0) = +--+ =0. The trivial solution 7, =2,=:- 
= n= 0 also has the foregoing initial conditions. In view of the unique- 
ness of the solutions of (1) under given initial conditions, it follows that the 
z; of (23) are identically zero in tf. 

Take any 7; =0 of (23) and differentiate the identity n —1 times with 


A, it js | 

1), 

lowing 
= 0, | 

3) the 
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respect to ¢. On setting ¢ 0 in each of these n identities, one obtains the 

m equations 

087 08’j10. 18j1¢3;°” > 1851404 + 0, 

(24) 2877 + 218’ 5162 + + ndjinCn™ 


(nm — 1) (m — 2) (m — 1) + AW" 


which are linear and homogeneous in c,, - - +, ¢n®. Since not all of the 
cx are zero, it follows that the determinant I of the coefficients must satisfy 
the condition 


0, 0, i, 1 

0, 1, Ai, Ae 
851875185 Sjin 2, 2A Ar’; 


which can be written in the convenient form 
(26) T= 5318 518518514 = 0. 


That the determinant I, is different from zero is easily shown, for, let A(A) 
be the determinant obtained by replacing the elements of the first column 
of T, by 1, A, A*, - - -, A** respectively. The equation A(A) —0 of degree 
m—1 in has Ai, Ai, As, An aS its n—1 roots. Consequently 


d?A(X) 
dd? 
for otherwise X =A, would be a triple root of A(A) =0. Since IT, ~0 at 
least one of the factors in (26), multiplying T,, must be zero. On omitting 
the corresponding column and considering the first n —1 equations in (24), 
one proves similarly that a second of the factors of (26) must be zero. By 
repeating this procedure one ultimately shows that 


=T, ~0, 


This result contradicts the fact that in each of the n solutions of (1) there 
is at least one non-zero coefficient. Therefore the assumption that A(t) =0 
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the is false, and hence the n solutions of Case II, Section 3, form a fundamental 
set of solutions of equations (1). 


0, 5. Derivation of a General Form of Solution for System (1). It will 
= 0), now be shown that when D(A,) and its minors have certain properties, the 
=i system of differential equations 


©) <9 has a solution of the form 
h (2) = + +4 + A,r) 
(t—1,:--, 
fy 
where the A,‘ are constants and , is a root of the characteristic equation 
D(A) =0. On substituting (2) in (1) and equating coefficients of t#,-- -, 
t° one obtains as sufficient conditions for the existence of a solution of the form 
(2) the w + 1 sets of equations 
2 
n 
(3) Ar) +S =0 (i= 1,---, ki), 
k=1 
y-1 n 
(4) (ais —,)Ai > = wAy™ (t= k t), 
k=1 
(5) (aii —Ar)AD™ + DS = (o—h + 2) 
k=1 


n 
k=1 


The solutions of these w + 1 sets of equations will be expressed in terms 
of the cofactors of D(A,) and their A-derivatives, and consequently it will be 
necessary to consider in detail the determinant 


Ay1—A, A102, A1,n-j» A1,n-j+15 Ayn 
Qe, n-js Qe,n-j+19 Gon 

D (A) == |Qn-j,1,  Un-j,29 An_j,n-j—A, |An_j,n-j+19 


— 
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where the jth and (7 —1)st principal minors, of orders n —j and n—j-+] 
respectively, have been set off with lines because of the important role that 
they play in the subsequent theory. On referring to D(A) one readily verifieg 


the identities 


(aii —A)8ij(A) + (2) (t=—1,---,n—j; 


(7) 


where the 5;(A) (t=1,-°-+, ~—j-+1) are the jth minor cofactors, re- 
spectively (of order n— 7), of the elements (1 =1,°° +, n—j +1) 
in the upper left (7 —1)st minor of D(A) ; and the 8;,;.,(A) ((=n—j+1, 

- +, m) are the upper left (7 ——1)st minors of D(A) obtained by replacing 
the (n j +1)st row by the (n—j-+1)st, (n—j-+ 2)nd,- - -, nth rows, 
respectively. 

For the foregoing particular choice in the position of the (7 —1)st minor 
cofactors of D(A) occurring in (7) the &; (k—=n—j+2,---, nm) are 
identically zero, by definition, and have been written for the purpose of 
securing generality in notation and logic, and also for the purpose of making 
the coefficients in (7) agree exactly, in form, with the coefficients in the left 
members of equations (3),---, (6). This specialization in the position, in 
D(x), of the (7—1)st minors in question has been employed solely for the 
purpose of attaining simplicity in presentation; all subsequent work, how- 
ever, is entirely general and applies to any minor-cofactors of D(A) satisfying 
certain conditions which will now be specified. 

Suppose (a), that (A—A,)* is a factor of every (7 —1)st minor-cofac- 
tor (of order n —j +1) of D(A), while at least one (7 — 1)st minor-cofactor 
does not contain (A—A,)**!; and (b), that (A—A,)* is a factor of every 
jth minor-cofactor of D(A), while at least one jth minor-cofactor does not 
contain (A—A,)**!, where obviously s—r=1. Let identities (7) be dif- 
ferentiated s— 1 times successively with respect to A, setting A= A., in each 
of the s sets of relations thus obtained. By hypothesis (a), the right mem- 
bers of the resulting s sets of relations vanish and, by (b), each term of the - 
left members of the first r sets of equations also vanishes. There remain, 
however, the following s—vr non-trivial sets of relations: 


(8) — Ai) + D = = 0 (t—1,---, 23 
k=1 
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(9) (ais — Ar) + (r + 
k=1 


k=1 
» 0; ky ‘ 
+1, It is to be remembered that at least 7 —1 of the 8,;,-- -, nj are identically 
cing zero in A, the remaining n—j-+1 of the 8;; being the jth minor-cofactors 
ows, of the elements of some row of a (j —1)st minor-cofactor of D(A,). Fur- 
thermore, the choice of the (j—1)st and jth minor-cofactors occurring in 
inor (7) is restricted only by conditions (a) and (b), where it is to be understood 
are that at least one of the 8;,;.(A) does not contain (A—A,)**1, and at least 
of one of-the 8;(A) (k =1,-- +, does not contain (A—A,)**1. 
ing If now the arbitrary » in (2) be chosen so that ow—s—r—l, it is 
eft observed that the w + 1 sets of relations (8),- - -, (11) have the same coeffi- 
in cients in their left members as the w + 1 sets of equations (3),---, (6). 
the The A,™,- +--+, An™ in the left members of (5) occur in the right members 
W- of the (h-+1)st set; corresponding relations hold in equations (8), °° -, 
ng (11). One thus obtains, by comparing the sufficient conditions (3),- - -, (6) 
with the equations (8), ---, (11), the following solutions of equations 
(3),°- +, (6): 
or 
(13) Ay? = C4; = Oo 
(h) wo(w—1) (w—h+2) (rth-1) 
(w—1) (o—h+2) coun 
(14) 


h+2 
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) wo! 
(w—1) ! 
C2; 2) 
2! 2) rt 


On substituting the values of the A,‘ given by (12), ---, (15) in 
(2), » particular solutions of system (1) are found to be 


= Art (r) 1 w(w—1) .(r+2) 
Ly = Cyje* to + (r-+1) (r-+2) 
wo! (8-1) 
1 (w—1)! 
J + [8x3 + r+1 8x3 t + + (r+1) (s—2) 
(16) 
2 
Aut (r) 42 (r+1) (r+2) 
[8x3 t + t+ (r-+1) (r-2) ] 
+ Caje™* ] (bel, n). 
After replacing the by (s—i)ey 
(t=1,+ + +,w), the solutions (16) can be written in the form 
84,5 -2) ¢ w-1) fw-1 | 


s—o—l1) ! 


The w particular solutions of (1) given by (17) can now be written in 
the convenient form 


wtl 


(18) caje* > 


(s—t) 
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Thus it has been shown that if, for any root A=A,, of D(A) =0 the 
(j—1)st and jth minor cofactors of D(A) have properties (a) and (b), 

then o =s—r—1 solutions of system (1) given by (17) or (18) can be 

written at once in terms of the jth minor cofactors of D(A,) and their A-deriv- 

atives. It remains to be shown that for every multiple root of D(A) =0 

the corresponding solutions of system (1), not of the type a —= Axe™* (k —=1, 
++, m), can always be obtained directly from form (17) or (18). 


6. Solutions of System (1) Corresponding to a Multiple Root of 
D(A) =0. Suppose (A) that A= A, is a root of D(A) —0 of order 
and (B) that all 1st, 2nd, ---, (q—1)st minors of D(A,) are zero, 
(¢q So), while at least one qth minor of D(A,) is different from zero. In 
consequence of (A) and (B), q purely exponential solutions of (1) are 
readily obtained in terms of qth minor cofactors of D(A,). For, consider 
again the system of differential equations 


(1) = HiT, Aintn (t—1,---,n), 


for which, if possible, solutions of the form 


(2) xi = AjeMt 


are to be found. On substituting (2) in (1) there results as sufficient con- 
ditions that (2) shall be a solution of (1) 


In view of hypothesis (A) and (B) the general solution of equations (3) 
can be written in the form 


(4) Aj = + C28ig2 + °° + Cgdigg (t= 1,°°°, 0), 


where the iq, are the gth minor cofactors of D(A,) obtained in solving (3) 
by determinants, and the c; are arbitrary constants. At least g—1 of the 
A; in (4) will have single-term values of the form Ai = cxdiqx. Substitute 
the A; given by (4) in (2) and one obtains the q purely exponential solutions 
of system (1) 


(5) = + + Codteqqe™* (k n). 


In consequence of assumption (B) all derivatives with respect to A, up 
to and including the (q¢— j —1)st, of all jth minor cofactors of D(A) vanish 


= 
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when A=A,, for j=1,-*-+,(qg—1). This follows from the fact that the 
A-derivative of every (j—1)st minor of D(A) is a sum of jth minors of 
D(A). Thus hypotheses (B) imply that (A—A,)%/ is a factor of every jth 
minor of D(A) [j—1, 2,-- +, (¢—1)]. If there exists at least one first- 
minor of D(A) which does not contain (A—A,)¢% as a factor, then the re- 
maining S) — q solutions of system (1), associated with A = 2,, can be written 
at once by (17) or (18) of Section 5, for the hypotheses of Section 5 are 
satisfied if 
S=8, r=—q—1, j—1. 


It may happen, however, that (A—A,)?, where p=q, is a common 
factor of all first minors of D(A) in which case there are fewer than S, —q 
solutions of the form (17), of Section 5, which are expressed in terms of 
the first minor cofactors and their A-derivatives. In fact, a similar special 
situation may exist for the jth minors of D(A) for 7—=2,---, (¢q—1). 
The assumptions with respect to this question will therefore be generalized 
so as to include every possible special case. 

Hence let it be assumed (C) that every jth minor of D(A) has the 
factor (A —A,)*’, while at least one jth minor does not contain (A —A,)S*1 
(j=1, 2,---, q), where obviously S; =q—j and 8;,—8S;=21. Hy- 
potheses (C) are identical in character with the assumptions underlying the 
derivation of forms (17) and (18), of Section 5, and consequently the results 
of the preceding article are immediately applicable. The solutions of equa- 
tions (1) thus obtained by successive application of forms (17) and (18) of 
Section 5, are, in general, expressed in terms of Ist, 2nd,---, gth minor 
cofactors of D(A,) and their A derivatives. In obtaining solutions in terms 
of jth minor-cofactors of D(A), (j=1, - +--+, q), one makes the following 
substitutions in (17) or (18), of Section 5: 


(6) s= = 8; r= §;. 


The total number of solutions, corresponding to X = A, which one obtains 
in this way, is clearly given by 


os (So — 8: —1)+(8: —8:—1) +: 
+ (Sq1—Sqg—1) = So — = So — 


since, by (B), Sg=0. There are also g solutions depending on the root 
A =4,, which are given by equations (5). Therefore there are S, solutions 


| 
| 
| 
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of System (1) associated with 4 =, which is a root of D(A) =0 of multi- 
plicity So. They are 


The solutions corresponding to every other root of D(A) =0 can be found 


similarly. 


%. Independence of the n Solutions of System (1). The independence 
of the n solutions of (1) will be established by proving that their determinant 
A(t) is distinct from zero. In order to attain complete generality in the 
following proof, assume that the characteristic equations D(A) =0 has B 


roots, Ax, * * Ag Of orders é2,° &, respectively, where clearly B Sn 
and + Let 
(2) = Ln = Pni (+=1, n), 


represent the n solutions of the system 

The determinant of the n solutions is 

(3) A(t) = | iz | = A(0) tome, 

The indirect method of proof will be employed. Assume 

(4) A(t) =0, 


which implies the existence of n constants ¢,,- + *, ¢n®, not all of which 
are zero, such that 


(5) 61 pir (0) din(0) n). 
From (5) it is evident that 

is a solution of system (1) having as initial conditions 

(7) 2,(0) = 22(0) =: -—2,(0) = 0. 


But 2, = 2, is also a solution of (1) having initial condi- 
t 
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tions (7). In view of the uniqueness of solutions (2) under given initia] 
conditions, (6) becomes 


(8) = ix +° bin 0 (+= 1, n). 
t 


Any 2; of (8) céa be written in the form 


where the f,i(t), 
are indicated by their second subscripts, while the coefficients occurring in 
these polynomials are either derivatives of cofactors in D(A;) or linear com- 


i=0 t 


-, fpi(t) are polynomials in ¢ whose respective degrees 


binations of such derivatives. The ny—1, nz—1, +--+, ms—1 are the 
highest degrees of the sets of polynomials fi:(t), foi(t),° fgi(t) re 
spectively. 


Set n, +m2.+----+n,—p and differentiate (9) p—1 times with 
respect to ¢, setting tO in each of the resulting p relations. One thus 


obtains the p linear homogeneous relations among the p constants ¢)‘),- + «, 


;(0) = 0, w’;(0) = (0) = 0. 


(10) 


Since (10) are consistent the determinant of the coefficients of the c;™ is 
zeTO, OF 


(11) T= II Hit Ti (a1, de) == (), 


i=1 


where 7;; is the factor common to the elements of the ith column in I, and 
is either a single term of the form om (Ai) or a linear combination of 
terms of that kind. In accordance with the foregoing notation se (Ax) 
indicates the ajith A-derivative for XA, of a pith minor-cofactor in D(A). 


After removing the factors from there remains the determinant 
(Ag?-1) 


where for convenience in writing 


ij 
i 
th 
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. In order to show that I, is not zero, replace the elements of the first 
column of T, by 1, A, A, respectively, and designate the resulting 
determinant by A(A). The equation A(A) —0 is of degree p—1 and has 
the roots Ai, As, * *, As Of orders ny —1, no, Mp Te.pectively. Since 
m—l+n.+:::+n,—p—1, all of the roots of A(A) —0 are ac- 
counted for. It is thus evident that 
d™-1A 
Tees 
: i 9, 
in 
— and consequently some 7;; in (11), say yj1, is zero. On omitting the last 
the equation in (10) and suppressing the terms of the first column of I, since 
) Te: ni: = 0, it is seen by a repetition of the foregoing argument that a second ji, 
th say nj2, must also be zero. By this procedure, one ultimately proves that 
wi 
thus (13) Wi: Yn * ‘= njp = 0. 
é since 7; was any one of the variables z,,- - -, 2n, it is thus evident that under 
the assumption A(t) = 0 one has (13) for 7—=1,---,n. 
In Section 6 it was seen that for any XA, the corresponding solu- 
) is tions of system (1) in terms of jth minor cofactors of D(Ax) have zero values 
for at least 7 —1 of the variables 7,,--~-, @%m. In view of this property it 
follows that for each A; there is an x; for which at least one of the correspond- 
ing nj of (11) is a non-zero term of the form Orne (Ai). Hence (13) 
contradicts the established properties of the solutions of system (1) and 
nd therefore the assumption A(t) —0, on which the derivation of (13) was 
of based, is erroneous. This establishes the fact that the n solutions of equa- 
i) tions (1) given in Section 6 form a fundamental set. i 
): il 
nt 8. Systems of Non-Homogeneous Differential Equations. Consider the f 
system of non-homogeneous equations i 
l 
(1’) aintn +O, (t= 1,---, 2), 
2 
: which differ from the system of homogeneous equations i 
4 
by the 6;(¢) which are assumed to be analytic functions of ¢. Let 
(2) = Gri(t), +, = Gni(t) (t—=1,---, 
represent a fundamental set of solutions of (1). Change from the dependent | 
4 i 
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variables 2; to new dependent variables y; by the transformation 


Since the $i; satisfy equations (1), the result of substituting (3) in (1) is 


On solving (4) for the y:, one has 


r 


(5) = Pity 95, > Din 
A(0)e™ 
where = dij n). 


On expanding the determinants occurring in (5), one can write the y’; 
in the form 


where the Pj;A(0) is the cofactor of pj; in the determinant | p;;|. Since t 
is the only variable occurring in the right members of (6), the y; are obtained 
by quadratures, and hence 


(7) > 6; Pye *tdt + (4 =1,---,n). 
j=1 


On substituting (7) in (3), one has as the general solution of equations (1) 


= =1 


A brief analysis of the nature of the solutions (8) will be given for each 
of-the following important cases: 


Case I. The characteristic equation D(A) =0 has no multiple root 
and (a), the 6; are polynomials in ¢; (b) the 6; are linear functions of Sin at 
and Cos at; and (c) the 6; are linear in Sin at and Cosat and polynomials in t. 


He 
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tiple roots which give rise to solutions involving polynomials in ¢, and the 6; 
have the properties (a), (b), and (c)_as in Case I. 


Case I(a). For this case the Pj, are constants and therefore the general 
solution (8) becomes : 


n n 
4,k=1 = 


After the indicated integrations are performed (9) can be written in the form 


where the y; are polynomials in ¢ whose degrees do not exceed the highest 
degree of ¢ in the 4;. 
Case I(b). The quadratures occurring in (8) are now of the form 


Sin at 
-r t 
at )e 


If a? ++ 42, 40 then on performing the integrations indi- 
cated in (8) one obtains as the general solution of (1) 


(11) Ry (Binet, Cotat) + Sosy (it, 
j=1 


where the R,; are linear in Sin at and Cos at. But if a? + A2, 0, then for 
Me = + aV— 1, 


Mat — 2atV-1 
J (Gos ( 1 [ ] 


For 4%—=—aV—1, the corresponding value of the preceding integral is 
obtained by replacing a by —a. The solution of (1’) now takes the form 


2, = R,(Sin af, Cos at) + hit (Cos af + i Sin at) + 
j=1 
(1=1, n), 


where the R; are of the same form as in (11) and the hy are constants. 
Case I(c). For this case the quadratures of (8) are of the type 


Sin at \,,, “hat 


CasE II. The characteristics equation D(A) —0 has one or more mul- 
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After performing these integrations the general solution of (1’) is 


(13) = H;(t, Sin at, Cos at) + > n), 


where the Sin at and Cos at enter linearly but have coefficients that are either 
constants or polynomials in ¢ of the same degree as the 6;. The special case 
for which 4» + a\V—1 makes no change in the form of solution (18) 
except that the degree of the polynomials in ¢ is increased by unity. 


CasE II. For each of the cases (a), (b), and (c), the Pjx of (8) are 
polynomials in ¢. In II(a) the integrals are of the same type as those in 
I(a) and consequently the solution of (1’) for this case is of form (10) 
where, in general, the degrees of the yi; exceed the maximum degree of the 6. 
Cases II(b) and II(c) lead to quadratures of the same type as those of I(c) 
and hence the solution of (1’) in either of these two cases can be written in 
the form (13). 
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On Generalizations of the Bernoullian Functions 
and Numbers. 
By E. T. BEtt. 


1. Both the numbers and the functions of Bernoulli can be generalized 
in an infinity of ways. Many of these generalizations have either an interest 
of their own or an importance in analysis. The methods by which general- 
izations are attainable can conveniently be segregated into two classes. Into 
the first fall all methods dependent upon modifications of the difference 
equations defining the Bernoullian functions. The recent generalizations of 
Norlund * belong to the first class. As we shall not again refer to this class 
it is unnecessary to describe it further. Into the second fall all methods 
dependent upon modifications of the symbolic equations of definition of the 
Bernoullian functions. The symbolic calculus implied is that of Blissard 
(and independently also of Lucas), in which the symbol A is the umbra, or 
representative, of the class of numbers or functions An (n 0, 1, 2,--°-), 
and A, is written A” until after the completion of all ordinary algebraic pro- 
cesses, when the exponents of umbrae are degraded to suffixes. During the 
intermediate steps powers A”, B”,- - - of umbrae A, B,- - - are treated as 
ordinary algebraic quantities, or ordinaries. For example, the first derivative 
of A” with respect to A is nA", umbrally, and this is nAn_, finally. We 
shall illustrate the second class of generalizations by an interesting example, 
which itself generalizes an extension of the Bernoullian functions and num- 
bers due to Krause and Appell. 


2. As the substance of Blissard’s long neglected papers is now well 
known, it will be sufficient here to recall a few essentials of his calculus as 
we shall use it. 

A function f(2™, - +--+, in which each (¢=-1,---, k) 
is replaced by the umbra A is evaluated as follows: the x‘ are replaced by k 
different letters A, B, +--+, C, each of which denotes an umbra, and after 
the degradation of exponents each of A, B,--~-, C is replaced by A. For 
example, 


(A + Aa)? = (A+ Ba)! = + A°Bta! = A, Bo + AoBya, 


*“ Mémoire sur les polynomes de Bernoulli,” Acta Mathematica, 43 (1922), pp. 
121-196. 
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if a is an ordinary, and the final value is therefore A,A,(1-+ a). Note 
particularly that zeroth powers of umbrae must always be indicated explicitly, 
In particular, A, B,- - -, C being umbfae, a, b, -, ¢ ordinaries, 

(aA +dB+:: 6C)° = Co; 
and the umbral multinomial theorem is, for n= 0 an integer (this theorem 
can be extended to negative and fractional exponents), 


n! 


the > extending to all sets (a, B,- - +, y) of integers a, B,- - -, y such that 

atBt:::+y=n. 
The umbral exponential theorem is, 6 being a parametric ordinary, A an 
umbra, 

AM /n! = And*/n! (0!—1). 

n=0 
This can be considered as an identity of definition. It follows from this 
definition that we do not have, for A an arbitrary umbra, e-4? equal to the 
reciprocal of e49, The value of e~4? is obtained from that of e49 by changing 
the sign of the parameter 6 in e49: 


(—1)"4,6"/n!. 
The one exception in which e~49 is the reciprocal of e49 is the trivial case in 
which A, =1 (n=0, 1, 2,°--). 

Infinite series, such as those just discussed, are used merely as generating 
functions, or generators. Thus e4® is the generator of An/n!, this being the 
coefficient of 6". The question of convergence of generators can be disre- 
garded. If this is not obvious, it has been proved in a former paper.* In 
particular we can equate coefficients of like powers of 6 in identities between 
generators derived by means of the umbral calculus. 

It is desirable to have no umbral derivatives in the final statements of 
results. Such derivatives can always be eliminated by repeated applications 
of the immediately verifiable identity 


— 


in which @ is an ordinary and A an umbra. 


* Transactions of the American Mathematical Society, Vol. 25, pp. 135-154. 
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Addition of umbrae v, p,- - -, ¢ is used here only in the following sense. 
If wis p,---, this is written u—v-+ p+ 
- + +-¢, and means that for all integers n = 0, 


the last of which is to be expanded by the umbral multinomial theorem. 
Thus if -+ p, then = + Vopr, U2 = + 20191 + Vopo, ete. 


3. Henceforth each of the k letters A, B,-- +, C represents an umbra 
of the Bernoulli numbers. 


Do = 1, D,=— D,= D,=— 1/30, - (n> 0), 


and each of the & letters X, Y,- - +, Z represents an umbra of the common 
Bernoullian functions with the respective ordinary arguments x, y,- °°, 2, 
so that 


1 

a r 

(W, w)==(X,z), (Y,y),°°°, (2,2). 

It is well known (or follows as the simplest case, k = 1, of the numbers 
© and functions A later discussed in this paper) that the symbolic equations 
of definition for D, W are 

0(1 — 
6 e9(1 — 


The arguments w—wz, y,* ~*~, 2 are not assumed to be integers unless so 
stated. 

We note first that Dn, Wn are functions of the integer variable n. A 
generalization which does not seem to have been considered is that in which 
D, W are functions of p integer variables &, yn, - - - ,€, 


W - 


that is, these D, W are members of p-fold sequences. A discussion of the case 
p=1 is a necessary preliminary for the case p>1. We shall therefore 
attend exclusively to p—1, reserving the similarly treated case p >1 for 
another occasion. 

There is a second possible generalization, applicable to p= 1, which we 
shall pass over here. We have postulated that the arguments z, y,---, z 
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shall be ordinaries. The case in which 2, y,---, z are umbral, so that 
x" = Xn, etc., is treated with equal facility by means of the theory developed 
in a former paper.* 

A third prolific generalization, immediately applicable to all that pre- 
cedes, consists in removing the restriction that F’, introduced presently, shall be 
algebraic in its arguments, a, b, 2’, y’,- + -, 2’. It is now sufficiently evident 
that the number of generalizations which can be systematically developed is 
infinite. The possibilities thus suggested will be more readily appreciated 
after the short treatment of the special case to which this paper is devoted. 


4, A power series in can be written in the umbral form 


where u,—n!tn. Let F=F(a, b, x’, y’,- ++, 2) be a rational algebraic 
function of its ordinary arguments a, b, z’, y’,- +, 2’, and in F’ make the 
substitution 

(a, b, 2’, +, 2’) = (6, €9, -, 


so that F becomes a function of 6, say f(@). Write the power series expan- 
sion of f(@) in the form e“@. Then w is an umbra of a system of polynomials 
in the variables z, y,- -, z, and = Un (2, y,* 2) is the nth 
such polynomial. For the special values (z, y,- 2) == (2, * 
where 2, Yo, * * *, % are given constants, the polynomials wu become a system 
of numbers 8. We shall say that the numbers £ are associated with the poly- 
nomials u, and that they have the characteristic determined by (20, Yo, 
* ++, 2%). Evidently u is uniquely known when F is given, and likewise for 
8B when in addition the characteristic is assigned. 

For simplicity let us restrict the discussion to the case in which the 
coefficients of F are rational, and the domain (of reducibility) is the rational 
domain. Then the resolution of the numerator of F' into irreducible poly- 
nomials in a, b, x’, y’,- + -, 2 is unique, and likewise for the denominator. 
Suppose common factors to have been removed. This unique form of F is 
expanded in the form e”?, and w is now called a primary umbra. From the 
definitions, primary umbrae arise only from irreducible rational algebraic 
functions F. The numbers associated with w are called primary. Non- 
primaries, either polynomials or their respective associated numbers, are called 
derived. 

A further generalization, easily attainable and readily developed by the 
symbolic method, at once suggests itself. In the foregoing we have replaced 
x’, ++, 2 by e%, respectively. But we might also replace 


* Mathematische Zeitschrift, Bd. 19, S. 35. 
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a’, 2 by eP@, eR), where P, Q,- +, are single valued 
functions. The case where P, Q,---, & are algebraic is of special interest, 
because it leads to an immediate generalization of the Bessel coefficients. 

It is now evident that we can set up a simple isomorphism between the 
multiplicative properties of rational algebraic functions in & + 2 variables 


a, b, x, y,* °°, 2 and the additive properties of umbrae of polynomials in 
k variables x, y,- * -, 2, and similarly for umbrae of the respective associated 
numbers. 


Suppose first that F is irreducible. Let e“# be the umbral exponential 
obtained from F as above outlined. Then from F", n >0 an integer, is 
derived e"-“9, where n.u is a convenient symbol for the sumu+v+-:--+¢ 
of n umbrae, each of which 1s to be replaced by wu after the degradation of 
exponents. Thus 


(3.u)? = (Ut + £)? = + Upvato + Uovote 
+ 2U1Vilo + Volz + 


If F is reducible, say F — G*H®- - - KY, let the respective umbral ex- 
ponentials of F, G, H,-- K be e%9, - where u, v, t,---, p 
are umbrae. Then clearly since for wu’, v’ arbitrary umbrae, e¥’%e”? = 
eu’ +08, we have 


With reference to umbrae, it follows (cf. § 6), that the rational algebraic 
functions giving rise to exponentials involving those umbrae play a part 
analogous to that of antilogarithms in ordinary multiplication. There is also 
a complete correspondence for division, but we shall not discuss this here. 

From the above sketch we see that the algebra of derived umbrae can be 
inferred from that of the primaries. The common Bernoullian functions 
A defined in (1) are primaries, since in this case F is F'(e%, e9”), and 
F(a, b, a, +, 2) =F(b, 2’) =b(1—2’)/(1—b), which is irreduci- 
ble in the sense defined. 


5. To illustrate the foregoing we take the immediate extension suggested 
by F(b, 2’), of §4 and multiply together & such rational functions. The & 
letters a, b,- - -, ¢ denote ordinaries; x, y,- - -, 2 are k independent variables, 
and as in § 3 each of the & umbrae A, B,: - -, C is an umbra of the Bernoulli 
numbers. 
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The Functions A, 


are defined by 


(2) eA@— 


With A it is evidently convenient to introduce the following, suggested 
by the right of (2): the functions of ¥, 


defined by 
(3) eve — (1— e276) (1 — ev?) (1 — e°9) ; 
the numbers Q, 

o*=0, 0), 


whose definition is 
ao bé 


the numbers 9, 
defined by 
(5) (1 — (1— - - (1— e), 
and the functions 9, 
= =F, ++, 2), 


whose definition is 


ax6 


Each of An, Yn, Qn, On, Bn, is of rank n, order k (=the number of 
independent ordinary variables z, y,- * -, z, or the number of ordinary para- 
meters a, b,: - +, c), and the one-rowed matrix (a, b,- --,c) is called the 
characteristic in each case; 7, y,* - *, 2 are the arguments of the functions. 

For the characteristics 


An is equal to Xn, Yn, - + *, Zn respectively, so that the common Bernoullian 
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functions are included in the special case / 1. For the characteristics i 
(a, 6b, 0,- 90), Cu. | 


we have 44(k—1) functions of order 2, of which a typical one of rank n is 
i 


! q 

so that the functions of Krause and Appell * are given by the special case 


= 2 with the further specialization that xz, y be integers =0. If 2, y, 
-, 2 are integers, the explicit form of A, is obtained at once from (2) 
rewritten as 


7=1 g=1 
so that | 
y 2 
The like is not true, of course, if x, y,- - -, z are not integers > 0. i 


From (4), (6) we have le 


and from (3), (5), | 
(8) b, (1, 1) — @,, 


For the respective characteristics 


becomes a"An, b"Bn,- c"Cn, where as always each of A, B,--:-, C is 
an umbra of the Bernoulli numbers. 


6. The relations between the numbers and functions defined in § 5 for 
various ranks are obtained by the following entirely general symbolic method. 


Let 6 denote a parametric ordinary, and 8, a, 8, - - -, y umbrae such that 
(9) by 

or what is the same by the definition of umbral addition, 
(10) (a+ (n=O, 1,2, 


* Ber. Ges. Leipzig, 55 (1903), pp. 39-62; Archiv Math. Phys. (3) 4 (1902), 
pp. 292-5. 
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Let the cn (n=0, 1, 2,-- -) denote ordinaries, and write 


(11) f(0) = 


where the upper limit r of the summation may be finite or infinity. If the 
latter, we shall assume for brevity that all infinite series occurring henceforth 
are absolutely convergent for some | 6|40. This assumption is not a neces- 
sity.* From (11) we have by means of (10) 


n 
and therefore 


(12) f(6+8) 


that is, from (9) we infer (12), where f(6) is arbitrary except that is pos- 
sesses a power series expansion (which may terminate) in 0. 

To obtain umbral identities such as (9), we observe that, if P(@), Q(@), 
- ++, 7 (0) are ordinary functions of 6 whose power series expansions in 6 
are given by the umbral equivalents ¢%, - -, eV? respectively, so that 


we have the symbolic equations 
P(8) =e, Q(6) = 0F,---, = 6%, 


and if moreover we have the ordinary identity 


(13) P(8) =Q(8) T(8), 
then 
(14) — +8, 


and hence we have (9). Hence the deduction of additive umbral identities 
such as (9) is reduced to inat of multiplicative ordinary identities. 
Consider next the frequently occurring modification of (13) : 


(15) KO'P(9) =Q(8) (9), 


where K is an ordinary constant (independent of @) and r an integer > 0. 
Then, instead of (14) we have 


— eat b+. -- 


*In the paper cited, I have shown how the contrary case is to be interpreted. 
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which can be written symbolically, 


4900 

08" 
and hence 

ar 

whence, as in deriving (12), and eliminating the umbral derivative by 
means of 

or % n 

(6+8)", 


we get, as the relation implied by (16), 
(16) K 7048) 


Hence by means of the simplest multiplicative identities between the 
ordinary functions defining umbrae we pass immediately to the generalized 
recurrences (12), (16). The process can evidently be applied as above to 
any functions $n, an, Bn, * yn (n=O, 1, 2,- which are coefficients 
in power series expansions.* ‘The coefficients themselves may be, as above, 
functions of several variables. 


?. By means of the processes in § 6 the relations between the functions 
defined in § 5 can be obtained with great ease. Change @ into —@ in (2) 
and compare the result with (2). Then, forn—0, 1, 2,---, 

[a(z+1) + B(y+1) =a, 
and therefore 
from which, by (12), it follows that 
(17) f(@—A+a(a+1) =f(6+A), 
A—a(x+1) —B(y +1) —e(2 +1)) =f(6—A). 


The special case in which A, degenerates to the ordinary Bernoullian function 
X,, is obtained from (17) by taking a = 1 and deleting the terms in b,-:-, ¢: 


* Tf, for example, R(6) = yp,6" defines p,, we first replace p, by n!p’,, so that 
> n n n n 
R(0) = yp’, n/n! =e 
and in the final results write p’, =P,,/n! 


i] 
| 

ta 
i 

| 
| 

| 

if 


286 BELu: On Generalizations of the Bernoullian Functions and Numbers. 


(18) 


the upper or the lower signs, as always, being taken throughout. For the 
functions of Krause and Appell only the a, b terms are retained from those in 
a, b,- ++, c, and further, x, y are restricted to be integers. Likewise in all 


subsequent formulas concerning A or the functions ¥, ®. 
In exactly the same way (3) gives 


(19) f(@+v) = (—1)* 
and (6) gives 


Proceeding similarly with (4), (5) we get for the numbers Q, ® the 
following, 


(21) 0), 


which could have been written down from (20) by taking in the latter 
y,° °°, %) == (1, 1); and 


(22) f(@+o) = (—1)* 
which follows also from (19). 


Corresponding to the relations between the common Bernoullian func- 
tions X, and the Bernoulli numbers B, for various ranks n, we have here the 
like for An, Qn, and the latter include the former as special cases. For (2), 
(4) with 6 replaced by — 6, (3) and (6) give immediately 

(23) A) = ¥ =O). 

Further general relations for A are obtained from (2), (5), 
(24) f(@+A+O—a—b—- -—c) =f(6+ 9); 
from (3), (5), and (2) with 6 replaced by — 8, 

(25) > Fez); 
from (2), (6), 

(26) +0), 
and from (2), (6), the latter with 6 replaced by — 6, 
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(27) f(@+A—®) 
From (4), (5) we have 
— (—1)*ab- = (—1)* ab: - cee, 

where «, = k!, en =0, nk, and therefore 
(28) f(@+2+40) = 

As an example of another type which, however, can be deduced from the 
foregoing by properly choosing f, let >. denote a summation with respect to 


all sets of 2 chosen from az, by,: - -, cz, a typical set of 2 being az, by, and 
similarly for >,, referring to all sets of r. Then from (2), and (4) with 0 


replaced by — 6, we have 


(—1)*ab---c — (1 — e278) (1 — (1 — e%8), 


whence, as in deducing (12), (16), we get 
ok 
(29) (—1)*ab---c age + A) 
0) —Dif(@=OQ+ az) + tart by)—:--. 
The similar theorem for Q is found from (4) rewritten in the form 
(—1)* ab: =e (1 —e%)(1— ¢%)-- - (1— e%), 
1 @ 
+3. 
The power of the symbolic method is well seen from the derivation of 
general addition theorems. Let (a’, b’,:--), °°) be the charac- 


teristic and arguments of Z’, (a, b”, ---), +) those of 
where Z’, Z” are not necessarily of the same order (same number of argu- 
ments where Z denotes any one of A, ¥, Then from 


(2), (3), (6) it follows that 


where Z is the function whose characteristic and arguments are (a’, b’,: - -, 
+ +), (a, y’,° +, 2”, y’’, respectively. Hence 
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(31) f(@+2) =f(6+2’+ 2”). 

If in the foregoing we take 2’ = y’ =- =y’ =: -=1, when Z is 
either ® or W we get 

(32) f(@+T) 


where 7’ is either Q or ®. 

As a last illustration of the power of this method let us find the partial 
derivatives of A, ¥, ® with respect to any of their arguments, say x. If Z, 
is a function of x, and Z’, its derivative with respect to 2, it is readily seen 
that the first of the following implies the second, 


0 
F(x, 0) = 0%, = F(x, 0) = 


To apply this to A, ¥, ®, take the logarithmic derivatives with respect to 2 of 
(2), (3), (5) and proceed as before, getting 


(33) (A+ + az)*, 
= (W+ aAz + az)", 
(35, = (@-+ aAz-+ az)", 


in which, we recall, A is an umbra of the Bernoulli numbers. Applied to 
(1), this gives for the common Bernoullian functions X, 


whence, integrating, we obtain 
1 
\ n+1 
(37) 1 (x A) 


it being readily seen that the constant of integration is zero. 

It may be mentioned in conclusion that the application of the symbolic 
method to Taylor’s theorem for functions of & variables, together with the 
preceding relations involving functions f of the several umbrae, leads at once 
to interesting generalizations of several well known expansions due to Euler, 
Boole and others in the calculus of finite differences. 
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The Construction of Certain Periodic Orbits of 
the Three Body Problem. 


By H. E. BucHAnan. 


In a recent paper * the author proved the existence of certain periodic 
orbits of three finite masses near the LaGrangian circular solutions. The 
present paper is for the purpose of constructing those solutions and discussing 
the geometrical character of the orbits. : 

The differential equations which define the variations from the circular 
orbits are 


— (1+ 8)*{(o? + 241)a, + 2Bixs} = (1 + 8)*e 
eU eU eU 

20(1 + — (1+ 8)?{(2B,2, + (o? + 243) a3} = (1+ 8)% 
dx 
t +8) — (1 (0? — Bays} = (1 + 8)%e 
d . 
+ 20(1 +8) — — (1 — As) ys — Bays} = (1 
+- 
d*z 
Sat (1 + Bits) = (148)%e 
Get (1 + 8)?(Bsz: + Asts) + (1 +8)’ 2x2; J 


* American Journal of Mathematics, Vol. XLV, Apr. 1923. 
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where the % denotes a sum of the form Ay,” + By,y; + Cy,” and all the 
constants have been defined in the paper referred to. 


Orbits wth Period 2r/v. We have demonstrated the existence of solu- 
tions of equations (1) which have the following properties: 


(a) Dy, Daj, (1—1, 3), 
j=0 j=0 
co 
= Daze’, Bi = Bie! 
j=0 j=0 


(b) Hach yij, 213, aij, Biz 18 uniquely determined, 

(c) 7% and 2,—c+y at 

(d) The series are convergent for | «| sufficiently small, 

(e) The whole solution and therefore each 2;;, yi;, 2ij, is periodic with the 
period 22/v. 


The coefficients of «°. The terms of the power series which are indepen- 
dent of « are determined by a set of linear homogeneous differential equations 
easily found from (1) by making e=—0. Their solutions satisfying condition 


ce are 


Lio = Kix 008 pit + Kiz cos wt + Kis cosh pot, 
Yio = Lis sin pil + Liz sin wt ote Liz sinh pot + Liat, 
Zio = Mi; sin ot + Miz sin vt. 


Where Ki; and L;; are uniquely determined in terms of L,; by the equations 


(A? — 2A,)Kyj— 2B,K3; 0, 
1; (A? — — 2A;) Ks; 2rAwLs; 0, 
Ky; (A? — o? A;) = 0, j= i. 2, 3. 


Let the initial conditions be 
Yyi= Bi, ty at t=0. 
The a;, Bi are related to the L,; by the equations. 
== + + 7 
a3 == + dg Lis + 
= Ly + oV —1 Lis polis + Lis 
Bs = Crp, V —1 Lis + —1 Lie + + 


(2) 


j 
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Where a,°, a3, B,° and 8; are the terms independent of ¢ in the power 


series which give a,, a3, 8; and QB; in terms of «. 
The periodicity condition (d) demands that 


Hence 


The M;; are related to the initial constants by the equations 


w? — A, 


oM,, + vM yo, 
3 


The periodicity condition (d) also requires that M,, = 0, hence 


| 
and 


c 


The terms of the series which are independent of ¢ are 
Lio = Yio = 0, 


210 == 


Gy sin vt, 


sin vt. 


aty 
30 = 
B, v 


The coefficients of «—These are determined by the following equations 
whose solutions must be subjected to the conditions (a), (b), (c), (d) 


and (e): 
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dy11 eU ) 
dt? — dt { + 2A1) + 2Bi%31} = 2"10 
dys1 
dé — 2w di {2Bst1 + (w + 2A3)Xs1} = 
eU 
2 
d*y dz @) 
Se —(w? — Ai) + Biysi = 0, 
a? = 
+ 20 + — (w? — As) yar = 0, 
d? 
dt? —_—-+ Aitu 28,{A, sin vt + B iM 3. sin vt}, 
de + + sin vt + 32 SIN vt}. 


The first four equations are independent of the last two. Their left 
members are the same as the left members of the equations which determined 
Lig and yio. The right members can be thrown in the form A, + A,.™ 
cos 2vt, Ago’ + Ase cos 2vt, 0, 0 where the A’s all carry the arbitrary 
¢; + y: as a quadratic factor. 

The solution of (53) consists of a complementary function plus a par- 
ticular integral. The complementary function has the wrong period as in the 
article just preceding and hence must vanish. 

The particular integrals are 


= Kio + Kio™ cos 2vt, = Lio + sin 2vt 


where Kio, Ki2™, Lio and Liz are determined so as to satisfy the dif- 
ferential equations. Applying the initial conditions we find 
== Kyo!) + Kig™, By t= 1, 3. 

In the last two equations of (4) the solution again consists of a complement- 
ary function plus a particular integral. But the right members have the 
same period as the complementary function. Therefore in order that the 
solution be periodic the coefficients of the right members must be zero. These 
coefficients are 


— 28:(AiMi2 + BiMs2), — 28,(BsMi2 + 
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2 
From the relation M3. —-t>“ 12 it is clear that the second factor 
1 


cannot vanish hence 6, 0. Since the right members vanish the particular 
integrals are zero. The initial conditions by (c) are = = 0 at t=0. 
Therefore the complementary function is zero and the solutions satisfying the 
conditions which have been imposed are 


om 0. 


The coefficients of «®. These coefficients are determined by the equations 


dy: 7 
{ (w? + + 2B,232} = 0, 
d 
— 20 + + 245) 200} = 0, | 
d dx 
— (a? — Ai) + Biyse = 0, 
dx. 2 | 
+ 2o + (wo? — As) ¥s2 = 0, (5) 
+ + ByZs2 
= — 28,(A,Mi2 + sin vt + sin vt + ¢,3 sin 3vt, 
dz 
282( + A3M 32) sin vt sin vt + sin 


The solution of the first four equations consists simply of the comple- 
mentary function. It has the wrong period and must vanish thus deter- 
mining as, and all of which are zero. Hence 


Vie = Yio = 0, 


are the solutions satisfying the conditions. 
In the last two equations the coefficients of sin vf must vanish giving 


— 28,(BsMi2 + AsMs2) + C31 = 0. 


Since Ms. = these give | 
1 
C1, By cs: 
2v andl 2v?(v? Mie 
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By direct computation one can verify that 


Hence the two values of 5, are the same. The right members of the last two 
equations of (5) will contain only sines of odd multiples of ¢. Therefore 
the complete solution is 


Zig = My; sin ot + Miz sin vt + K; sin 3vt, 


where the K;‘? have been determined by the usual method for finding the 
particular integral. By the initial conditions 


o — Ar v A, My. + 3vK, = y3. 
B, B, 


Hence setting M,, = 0 as the periodicity conditions require we get 
A, 
B 


1 


A direct determination of the K’s shows that this value of y,‘* is not zero, 


The Normal Form for the Z-Equations. In each step of the construction 
there will appear two determinations of 5;. In order to show that these two 
values are always the same we make a transformation of the form 


+ V2 + V3 + 

= + V2) + + V2), 

=wV —14,(v,— v2) + vV —1 d2(v3 —%). 


In the normal variables the Z-equations of motion become 


(1+ 8)oV —1 0, + (1+ + CRs}, 

=—(1 + 8)oV —1% + (1+ 8)e{— CLR, — C3R3}, (7) 

== (1 + + (1+ $)e{— + 

=—(1+ 8)vV —1 0, + (1+ 8)e{C,Ri — 
Where FR, and FR, are the right members of the original Z-equations. It is 
not necessary to change them into v’s. 

If we had made the previous part of the construction in the new variables 

it is clear that the coefficient of e”V-1 * is the same in the vs and v, equations. 
Hence the two ‘determinations of 8, are obviously the same. 


A; (2) 
Cy) 071. 
B, 
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The general term. Suppose that the 
have been determined and have the following properties: 


(1) Liok = Yiok = Zioke = = 0, 82x 


(2) Dioxy1 are cosines of even multiples of vt, the highest being 2(k +1). 
(3) Yioxsr are sines of even multiples of vt, the highest being 2(k +1). 


(4) Ziox are sines of odd multiples of vt, the highest being 2k + 1. 


It is possible to make an induction showing that these properties persist. 
This was done by Moulton in his Oscillating Satellite paper, for the case of 
an infinitesimal body. The details are so closely related to his work that it 
is not necessary to repeat it. 


The character of the orbits. In the generating solutions 7; = y; = 0, 
z,—= M;sinvt there are apparently two constants M, and M; but Mj=—= 
A, 

B, 
solutions show that c, + y, enters as a factor of each Zin, Yin and Zin but 
to a power one greater than «. Therefore, if we replace (c, + y:)e by ¢ 


M, so that there is but one independent constant. The forms of the 


we get 
oo co co 
(1 +n) rij, Yi= (er + 71) (+71) 
=0 j=0 
The coordinates in the orbits are given by 


j=0 j=0 j=0 


where the primes denote the rectangular coordinates in rotating space. The 
exact forms of 2’;, y’i, 2’; as far as the second power of e’ are 


Qvt 
2’; = My sin vt + J 
1+8 , 


The equations of the projections of the curves on the plane is found by 
eliminating ¢ between the first two equations. The result is 
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— Kip 4 4 { y's 
‘ 


which are the equations of ellipses whose centers are on the z-axis. 


The projections on the 2’, 2’ are the parabolas 


22/2 


The projections of the y’, 2’ plane are 


The orbits are therefore curves shaped like the figure 8. It can be proved 
from equations that the parabolas are opening towards each other and that 
the values of y; for small values of ¢ are negative. The orbits are therefore 
retrograde. ‘These projections show their general shape with reference to the 
rotating axes. The orbit of the middle body is always determined by the 


center of gravity equations. 


Orbits with the period 2r/w. The properties of the solutions which were 


established in the existence proof are: 


(a) each 2ij, yij, 2ij, 8; is uniquely determined. 
(b) = 2, = 0 and 2, + at 
(c) The series are convergent for0 << 


(d) The whole solution and therefore each 2;;, yij, zi; is periodic with 


period 27/w. 
Applying these properties we find the coefficient of e° to be 


Lio = Ki, cost, yio= Li sinot, = Mi, sin of 


since by the condition (d) the coefficients of cos pit, cosh pot, cos vt and of ¢ 


must vanish. That is, we have 


Li = Lig = Ly = Mig = 0 


These equations are no longer homogeneous as in case 1 for the initial 


conditions are 


The first three equations have a determinant different from zero. 


| 
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it is possible to solve for three of the initial constants in terms of the fourth, 
which was taken arbitrarly, and without loss of generality can be made equal ~ 
to zero. The equation M,,;— 0 gives y; uniquely in terms of y; which is so 
far undetermined. It will be necessary to determine y, along with 8 so that 
the solutions at the successive steps shall be periodic. 

The coefficients of « are determined by the differential equations. 


— — (w? + — ) 
= 281 { (w? io +- 2Bixxo} + + Ain cos 
28, {w? Ai) Yio — Biyro} + Bi.™ sin 2ut, (9) 
= — 28, (Aizio + Bizeo) — yi (Aizio + Bizeo™) 
+ Co + Cis sin (i—1, 1). 


The superscript notation here means that the part of zo going with the first 
power of « in y; is taken. In these six equations the first: terms of the right 
members have the same period as the solutions of the left members set equal 
to zero. In order that the solutions shall be periodic we must have 


=8, =0. 
Then the solutions become 


= Aio™ COs wi cos 2wt, 
= Bio? + sin ot + Bi 3in (10) 
Zi1 = Cin? + sin wt + sin 


where the Aij, Bis, Ci; are determined so as to satisfy the differential 
equations. 

In the computation of the coefficients of <? the equations which determine 
yi and 8, are no longer homogeneous as above. Further, there are six linear 
equations which these two quantities must satisfy in order that the solutions 
shall be periodic. In order to show that these six equations reduce to two 
independent ones it is convenient to transform to the normal variables for 
the remainder of the construction. 

The equations of motion become in the new variables.* 


* See equations (22) Am. Journal of Math., Apr. 1923, p. 101. 
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w= (1+8)6V —1 wt (148) Anel +<«( 


aU eu 
+) 


=— (14+8)oV —1u+ (1+ 8)e 
u's = (1+ 8) + (1+ d)e 
=— (1+ 8) + (1+ 


= (1+ )d ) + ) }; 
= Ur + (1+ 8)e 


+ (1+ 8)?(Aiz, + Bizs) = (1 + 8)%e 


+ (1+ 8)?(Bs21 + Aszs) = (1+ 8)%e 2x2; 
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+ e( ): 


+e( 


= — (1+8)d:V + (1+ 8)e 


u’s = (1+8)oV + (1+ | (11) 


{Bu( ) + ) + Bis( ) + Bul )}, 
{— Bu( )—Bi( ) —Bis( ) —Bul 
{Cu( ) +Cr2( ) +Cis( ) + )}; 


{—Cu( ) —Ci2( ) —Cis( ) — Cua )}, 


{D1:( ) + ) — ) — )}; 


where the parentheses denote the coefficients of Ai:, Aiz, Ais, Ais in the first 


equation. 
For the sake of clearness we will compute again the coefficients of e° and 


e. Let the initial conditions be 
+ yt, ci = 


Therefore the coefficients of «° are 
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Orbits of the Three Body Problem. 


Uso = (4 + a3) #, (a + ay) 
Zio = Min sin wt. Uio = 0, (*=1, 2, 5, 6, rf 8), 


By the orthogonality conditions a; = a, and therefore 


= a) cost, = + a) COS wt, 
Yio = 2(a + a) sin wt, Y3s0 = 2c, (a + asin ot, (12} 


Zio = Mio Sin ot. 


The coefficients of « are determined by a set of equations like (11) in 
which all the terms containing e to a higher power than the first are omitted 
and the values of 2i0, Yio, Zio are substituted from (12). It is easily seen 
that the periodicity conditions are satisfied by 6, = yi" —0. The solutions 
are 


= a, + sin 2wt + as cos 2vt, 
Uo, = —a,™ sin 2wt + a3 cos 2wt, 
= do + gin + bs cos 2ut, 
U4 = t 4+ — gin 2wt + cos Qt, 
gin 2wt + cos 2ut, (13) 
Us. = — sin 2wt + cos 
+ d,™ gin 2ut + ds cos 
Us: = — d,™ sin + d; cos 2wt, 
My” sin wt C43? sin wt. J 


The two arbitrary constants of integration, bo and 6” are determined by 
the orthogonality conditions so that bo’ —b® —6/2 where b is a known 
quantity. Therefore in the x, yi, 2; variables we have 


= Ato’? + cos ot + Aio™ cos 
4a = Bio Buy sin wt + sin 2wt, 
= Cio? + sin ot + sin 2ot. 


The character of the differential equations together with the initial con- 
ditions shows that 2; will always be a cosine series, y; a sine series, 2; a sine 
series, provided y,‘"5, can be determined so that the solutions are periodic. 
We will prove that there are just enough independent equations to determine 
yi™ and 8, uniquely. The equations defining the coefficients of e* are 
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= OV —1 tan + —1 + Ago™ + sin ot 

+ sin - -+ cos wt + cos 2wt 
Won = —1 tan + V —1 + + Asi ™ sin wt 


2” in = — — ByZgn — + — yi™ + ) 

+ + sin wt + sin +: - - 


where the equations omitted are already in such a form that the periodicity 
conditions are satisfied. The condition that the first equation of (14) shall 


have a periodic solution is 


since the coefficient of e*V- * must vanish. In the second equation the coeffi- 
cient of e’V-?* must vanish, which gives exactly the same equation as the 
above. This equation uniquely determines $n. 
In order that the last two equations shall give periodic solutions we must 
have 
28n(AiMio + BiMso) + (AiMio + go) = $12 


These two equations can be shown to be identical by reducing the z; equations 
to the normal form as in the last article. Since 8, is already known, y,‘" is 
uniquely determined by either of these equations. 

In the generating solutions there are two independent arbitrary numbers 
L,, M,. It is no restriction on their independence to put M;—CL,. Then 
in the solutions constructed above, Z, enters as a factor of each xij, yij, 24; 
but to a power one greater than «. Therefore in the expressions for the co- 
ordinates, LZ, can be absorbed into «. The C is left arbitrary. 

We have investigated the orbits with the period 27/w under general 
initial conditions, and under orthogonal initial conditions. Let us now con- 
sider the possibility of orbits with period 27/w such that the coordinates are 
in the following proportion 


x 2 


i 

| 

| 
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If we make this substitution in equations (2), they reduce to the two-body 
problem and are completely integrable.* 

When the above proportion is satisfied the bodies all move in ellipses and 
remain in a straight line. The plane of motion may have any inclination 
whatever, which accounts for the arbitrary C left in the construction above. 


Orbits with Period 2r/p;. We have already proved that these orbits lie 
wholly in the zy-plane. The method of constructing them is very similar to 
that of case 2 with the z-equations omitted. It is again necessary to trans- 
form to the normal form in order to show that the two determinations of 8, 
are identical. 

The geometrical character of the orbits is described by the following 
statements : 


(1) The orbits are symmetrical with respect to the z-axis geometrically 
and the time. 

(2) Omitting powers of e beyond the first, the orbits are ellipses with 
centers on the z-axis. 


Orbits in which p, and v are commensurable. These are of the same 
character as the orbits with period 27r/w. The same difficulties of construc- 
tion arise and may be handled in a similar fashion. The orbits are symmet- 
rical with respect to the z-axis geometrically and in the time. In case (2) 
we showed that the orbits lie wholly in one plane with an arbitrary inclina- 
tion but in this case there does not seem to be such a simple interpretation 
of the constants. The projections of the orbits on the zy-plane; neglecting 
powers of « beyond the first, are ellipses with centers on the z-axis. The 
period is the least common multiple of 27/p; and 27/v. 


* Lagrange, Collected Works, VI, pp. 230-270. 


Correction. 


On page 223 of the current volume of this JouRNAL, the first parametric 
equation for curves (9) and (10) should read 


= + at. 
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